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ABSTRACT 


é- 


It is known that the static and dynamic behavior of 
ieuids in Peres wmMediasdepends TO amma se measure Gn poreus = 
media geometry. In the past,.the abittry lo Characverize 
this geometry has been restricted to such average properties 
as porosity and permeability. However, in recent years 
attempts have’ been made to achieve a more precise charac- 
terization based upon the fact that porous media is 
statistically composed. 

In this thesis techniques from statistical ‘communi -: 
cation theory are adapted as possible methods for accomplish- 
ing this classification process. In simplest form, a 
euchotamous function is defined by passing a line throughea 
porous mediun, the function having one value when the' line 
is in solid matrix, and another value when the line passes 
through pore space. The function is then analyzed neAte 
(1) Classical Fourier series harmonic analysis, and (2) 
determination of the autocovariance estimate and power 
spectrum. Comparisons are made between several functions 
created from the same medium, and with functions created 
from other media. 

The results indicate that the autocovariance estimate *° 
and the power spectrum, as characterizing functions, oth 


ilps leat 





discriminate between different media. This success suggests 
many more paths of investigation, possibly leading to the 


complete characterization of porous media through statistical 


analysis. 
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CH APT Era 
INTRODUCTION 


In August, 1965, two American astronauts journeyed 
over two hundred miles into space in preparation for man's 
future journeys to the moon and thousands of miles beyond. 
As of August, 1965, man had journeyed only a few hundred 
feet beneath the surface of the earth upon which he lives. 
In this very interesting paradox, man is now capable of 
gathering information many miles into space, but has very 
little first-hand information about the sub-surface of the 
very earth from which he launches his space probes. If this 
seems incongrous, consider that while the benefits of space 
exploration are speculative at best, man derives the vast 
majority of his resources, his energy saurces, and even his 
food from this earth. However, there are many explanations 
for man's actions. Perhaps the foremost is the extreme’ 
difficulty encountered in any attempt to explore extensively 
beneath the earth's surface. Costs would certainly be pro- 
hibitive, even if the engineering difficulties could be | 
overcome. Another explanation might be that man has been 
able to theorize sufficiently about the earth's sub-surface 
areas from the limited information available so that he has 
not been willing or seen the necessity of risking great 


di 





expense, for what might be very little gain. A less prac- 
tical explanation, but in all probability ea ver Beedle. on 
.for the lack of extensive sub-surface exploration, would be 
that there is less glamour to most people about traveling to 
the center of the earth than in a journey to Venus or Mars 
and thus, less interest in making the attempt. 

The above discussion is not meant as an implication that 
man is not interested in the earth's sub _cuRece: To the 
Cenlrary, @eOkOpiste, petroleum engineers, anameuicnr. who 
derive their livelihood from the soil and beneath the soil, 
have performed extensive research in this area. From the 
exploration of caves and mines, from core samples drawn from 
wells, and in many oonen pate man has been able to learn 
much about the earth's constitution and composition. Asa 
mesult of this research, great improvements have been realized 
in the methods used to extract the natural resources from | 
the earth, not to mention the improved conservation procedures 
which will allow man to benefit from the earth's resources 
for many Venue to come. | 

As a part of this research, man has spent many years 
aeveloping intelligence concerning the physical properties 
of porous media.found in the strata of various kinds that 
Jon: pence earth's surface and the behavior of fluids 
within these porous media. Research efforts have been con- 


centrated on four properties which are considered fundamental 


and the basis for discussions regarding all other properties.~ 
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These four properties are porosity, a measure of the void 
Space in the media; permeability, a measure of the media's 
ability to permit flow; the fluid saturation, and the 
electrical conductivity of the media. We will discuss these 
properties at greater length in Chapter II. 

The basis for this report is: ‘derived i vom a papemsoy 
eee Scheidegger“ on "Statistical Geometry of Porous Media," 
published in 1961. As wiuaees out by scheidegger, although 
the usual geometrical quantities, such as porosity and 
permeability, give generally good descriptions, there is a 
need for geometrical description of greater detail. This 
need arises because of the extreme complexity of the subject 
matter. For example, in some cases, media with a given 
porosity will not allow as much flow as other media with less 
void space because of the lack of continuity in the pore 
spaces of the first sample. Media with identical permea- 
bilities to one fluid may have different permeabilities to 
another fluid. These are well-known phenomena in the area 
of petroleum engineering. As a result, it is sometimes very 
difficult to classify media of different types without 
relying on extensive laboratory tests, which are expensive 
and time consuming. Po wts viega) teen ens report, to take 
advantage of the inherent statistical nature of porous media 
to characterize the media by a set of numbers or a mathe- 
matical ae ene To clarify the above statement, Vigeains 


important to point out that porous media are consolidated in 





mM 
a statistical nature. An individual» pargiclespaewp oe ler. 
determined by its size, shape, density, electrical conduc- 
tivity, and various other characteristics. The total media 
Will have properties determined by the properties ol vuie 
Meaividual particles, and by the means in which the individual 
Marvicles' hayé united to form the wholes) This cConscligauien 
process will in turn depend upon the arrangement of the 
elementary particles, the relative nee Of Gilftecrent types 
of particles, and upon the relative number of large and. 

small particles. This dependency upon Gletoreunen [UNC los 
is obviously statistical in nature. 

While much time has been devoted to research on the 
various physical properties mentioned above, very little has 
been accomplished in the latter area of investigating the 
statistical nature of porous media. Scheideggar? and Flood 
have provided impetus in this area, and defined some of the 
Foals. It would be particularly desirable, for example, sega 
be able to look at two sees of porous media through a 
microscope and tell if they are identical in a statistical 
sense. It would be even more desirable to be able to make 
intelligent predictions about the other fundamental properties 
from the statistical description. 

The purpose of this report is to introduce and watcemps 
techniques which utilize statistical communication alee) 1, 


characterize porous media. In simplest form, the method will 


consist of defining a function created by passing a line over 





5 
or through a porous medium, and assigning a value of *1.when 
the line passes through pore space, and -] when ice passes 
through solid matter. This function will then be analyzed 
in a number of ways. One attempt will be to fit the function 
with a single Fourier series and determine whether the 
Fourier coefficients effectively characterize the medium. 
Brief mention will be made of the double Fourier series 
whereby a surface area, rather than a single line, might be 
analyzed. 

Finally, the pict tion will be looked at in terms of cer- 
tain statistical properties, derived from communication | 
theory. The idea of random functions will be introduced and 
elaracterizing tunctions such as the autocovariance TUunGGaen 
and the power spectrum will be determined. Comparisons will 
be made within the same medium, and against different media. 
From the results, it is hoped to ascertain whether a given 
porous medium can be described and effectively classified by 
its power spectrum and autocovariance function, to the: 
exclusion Or all others: 

If any of these three methods can be used to effectively 
Classify porous media, the obvious step to be taken later 
will be the Sonate of these A anonetee quantities with 
the fundamental properties, such as porosity and permea- 


bility, which are now used to characterize porous media. 





CHAPTER Ii 
HISTORY 


Until quite recently, man has investigated and classi- 
Brea porous Media largely according .tco fouges Undameiua 
properties from which it was felt all other important prop- 
erties could be derived. These four basic properties are 
permeability, porosity, fluid saturation, and electrical 
conductivity. sSome of the earliest work was done by Henry 
Darcy, a Frenchman. In 1856 Darcy, investigating the flow 
of water coe a sand filter, equated flow rate with the 
cross sectional area of the conduit, the length, and the 
pressure head Creaveay oy having the ends of the conduit at 
different lene He found that a Soporte enon constant 
was required to balance the equality for various filters, and 
that the constant was a function of the sand packing. This 
constant became known as the permeability of the media. The 
unit, of permeability is appropriately called the darcy, and 
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is defined in the petroleum” industry as: 

A porous medium has a permeability of one darcy 
when a single phase fluid of one centipoise viscosity 
that completely fills the voids of the medium will 
flow through it under conditions of viscous flow at a 
rate of one cubic centimeter per second per square 
centimeter cross sectional area under a pressure or 
equivalent hydraulic gradient of one atmosphere per 
centimeter. 





( 

As can be seen from the definition, Darcy's law may be 
extended to include other fluids as well as wate Silene, 
the permeability is a function only of the porous medium 
and not the particular fluid flowing. However, investigators 
such as she ernagpe reported permeability variations dis- 
covered when using gases as the flowing fluid. The permea- 
bility was found to vary with the gas pressure, and the 
variation was attributed to a phenomenon known as gas slip- 
page. A discussion of this phenomenon is presented in most 
standard texts, e.g., Amyx, et Bg! 

Prior to discussing the work of Carman and Kozeny in the 
prediction of the permeability of a porous media from its 
_known pore size aeteaoun tear we should discuss early work in 
the area of calculating porosity. The early investigations of 
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Blichter” and Fraser and Graton’ were concerned with the 
Peeadicvion Of porosity of packings of uniform Spheres. Since 
porous media are, in most cases, quite obviously not composed 
SueuUniform spheres, but of a myriad of particle sizes and 
shapes consolidated with every meceie ee angularity, the 

work of these men must be enlarged upon to deal correctly with 
real materials. Whereas their work was largely theoretical, 
and others continued to work with predictions of porosity based 
upon grain size distribution utilizing seive analysis, the 

more recent work in this area has been through empirical 


determinations of porosity, by laboratory measurement, using 


such devices as the Kobe porosimeter?® (mercury introduction), 





ora the Washburn-Bunting? porosimeter (air extracted by 
vacuum and measured). 

studies of pore size distribution have been carried out 
by several investigators including Carman, ?* using gas 
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adsorption techniques; Ishkin and Kazaner, passing two 


phase fluids through the media; Ritter and Erich,?¢ using © 
small angle x-ray scattering patterns, and Sy anemone 
with variety of techniques. Grain size distribution work 
has been accomplished mostly by simple sSeive analysis and 
there is apparently a need for more work in this area. In 
fact, Masson? indicates that after extensively researching 
this field, there is "a definite need for an inclusive 
mathematical function containing a few definite parameters 
that would define a given size distribution for most con- 
ditions approaching ideal grinding." 

One of the most interesting works done in the deter- 
mination of pore size distribution is attributed to Nuss and 
Whiting,2° who devised a method of impregnating core samples 
with plastic under high pressure. The plastic was then 
polymerized at a high temperature and the rock material 
removed by acid leaching. The remaining plastic reproduces 
the original DOTS  Ssrace. | 

Many investigators have completed correlations between 
the porosity and permeability of porous media. Carman, ~! 


Kozeny,—~ and Wyllie? developed predictive equations from 


these correlations, whereby one property might be estimated 





accurately from the other. In recent years in England, 


ae have refined the Kozeny Equation 


Pollis-Georgesand Childs 
to improve results and allow greater generality. Francher, 
Lewis and Barnes“ experimented with various porous matter 
pe debermine vwne relationship bétween grain suze ange 
conductance. In their work they accounted for both laminar 
and turbulent flow. A brief discussion of the results iL 
found in the text by Amyx .°* 

Nearly all rock samples brought to the surface contain 
amounts of liquid entrained in their pore space. The deter- 
mination of the amounts and types of liquid is a very 
important characteristic description of the porous media. 
Comprehensive descriptions of the methods used to accomplish 
this determination are covered by eicigti ya Hornene 
Schilthius,“° and others. The two most common techniques 
for determination of fluid content are by retorting, and by 
solvent prract ln and distillation. Fluid saturation may 
also be determined by using an indirect approach, that is, 
by measuring some other property and correlating the asa Bkiae 
with ee data to estimate the fluid entrained within the 
Sample. Use of electric logs and capillary pressure measure- 
‘ments are two such indirect: means of arriving at the fluid 
saturation. 

The fourth fundamental property of porous rock is the 


electrical conductivity. The rock, @s 10 Comes Tremeune 


ground, is comprised of solid fragments, void:space and some 





dee 
amount of entrained liquid. Generally, the solids and the 
Oil and gas eonponents filling the pore spaces are non-— 
@onductors of electricity. The wateryeit it contawisea.. 
solved salts, then, is the only conductor. The ability or 
non-ability of a rock to conduct electricity is thus an 
madication or the HE eMiegeteze} a SKS 101553 or continuous water filled 
pores. Wyllie and Speynetlens © created the firs. modelwon 
porous edie taking advantage of this phenomenon to charac-.- 
terize the media according to its electrical resistivity, 
and hence, derived such other properties as tortuosity and 


a and Wyllie and Garaner“® 


Porosity. Cornell and Katz 
added refinements - later models. 

From the four basic properties covered, many other 
properties can be derived. Researchers in many fields have 
developed various descriptive properties which characterize 
fede more explicitly for their particular use. For example, 
the geologist has expanded the characterization to include 
mineral content, while the petroleum engineer, because of his 
interest in flow properties, has developed such concepts as 
capillary pressure, relative permeability, and wettability. 

A description of these properties is beyond une scope of this 
report, bue the bibliography references the works in this 


‘ 


area of patel <2 everett. Flee Brown, + Rose,“ Slobod ,?- 
Fatt," and Rap@port and Leas .?” ° 
Man has, then, taken advantage of many physical properties 


to descriptively classify porous media. However, in 1961, 
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Pe. oChermaeeseer Andie a lime icmess Writing. nee ale irae ona 
Geophysical Research, point out that the usual geometrical 
properties, such as porosity and permeability, do not describe 
the media sufficiently for present requirements. Quoting 
Eneir article: | 
Geometrically, a per cucenedetin is det inedmoy 
giving the analytical equation of the surface which 


bounds its pore space. For any practical purpose, 
this is impossible to accomplish. Average geometrical 


@uantities, such aS porosity. ~~ | aiddearce ae 
surface... , are conceptually acceptable (and can 

be measured), but this is not true for the concept of 
"pore size." Since the pores form an extremely compli- 
cated system... , it is not easy to determine a 


pore size in a conceptually satisfactory fashion. 
Nevertheless, there is a great need for esta- 

blishing a better geometrical characterization of a 

porous medium than is afforded by the mere deter- 

mination of porosity and specific surface area. 
Scheidegger and Fara desired a means of characterizing porous 
media through analysis of the statistical make-up of the 
media. | 

These investigators introduced four possible methods of 
statistical analysis: (1) by use of a correlation function, 
(2) a spectral analysis in terms of harmonic functions, (3) a 
spectral analysis in terms of other orthogonal functions, 
and (4) a spectral analysis of a specially constructed éhancem 
teristic function of the media. 

In each of the four techniques, the writers developed 
statistical data for interpretation by drawing an arbitrary 


line through the porous medium. The line has a designated 


origin and data points are selected at equally spaced points 





ky 
allow tne line. At each of these points, a value of *1 1s 
assigned if the line is, at that point, passing through void 
or pore space, and a value of -l is assigned if the line is 
passing through solid. This series of *t1 and -1l values, 
taken from the origin, define an orthogonal function f(s), 
which is characteristic of the line drawn through the funct1 or 
and of the media itself. If the medium is homogeneous and 
meourOpic, andeift the sample length isvof sufficienussize. 
theoretically, the functions or lines are independent of 
direction, origin, or position. The function is treated as 
a random function, as defined by Moyal,?¢ and can be statis- 
tically analyzed in several ways. 
If the function is defined as above for void and solid 
spaces, the mean of the function may be derived: 
Perle (-1) (lee) ee) 
where P * porosity or pore spaces 


f er = |. | (2 2) 


In many bacon, scheidegger points out that it is anecue 
venient to work with a function having a mean other than zero, | 
so he rewrités (2) as a new function with a mean that is 
equal to zero. The new function becomes: 

7 F'(s) = f(s) - (2P-1) (2.3). 

In the first of the four vee aS which the authors 
suggest, the function is charac temiees by its aU LOC ORES aman 


PumCcit Tem. 





= . us 
R(A) lim ‘ae les f(stA)ds (2.4) 
S 


a ds 


Chapter III of this paper SOntalnome discussion of the auto- 
Serrelation function. 

The second technique referred to in the Scheidegger 
‘paper is to fit the function, f(s), with a Fourier series, 
and evaluate the Fourier coefficients, a, and DL. A suitable 
evaluation suggested is to use the equation: 


Dm 
a 


aye 2 1p (2.5) 
where a is not only a fair representation of the medium's 
geometry (statistically), but has the added advantage of being 
macependent of point of origin or phase angle. 

The third and fourth techniques suggested are mere 
refinements of the two discussed briefly above. Each of the 
ie War LeCnTGMes was offered by Scheidegger and Fara in an 
introductory manner, and there was very little in the way of 
tangible conclusions. However, these men opened an entirely 
new facet to the subject of characterizing porous media. 
Unfortunately, in correspondence between the author of the 
report in hand, and Dr. Scheidegger, the latter indicated that 
he had done no further work in this area. ee works ‘in 
the refinement of these and similar techniques are now left 
to this author and others. 

Perhaps brief mention should be included herein as to 


the history of statistical communication theory and time 
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series analysis, which Scheidegger has introduced into the 
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world of porous matter. Granger credits the astronomers 
as being the first to analyze time series, although he men- 
tions discussions of meteorological and geophysical series 
in the middle of the nineteenth century. Even in these early 
discussions, it was noted that series of data taken over > 
periods of time often follow patterns and contain certain 
trends, sometimes obvious and sometimes not so obvious. If 
the obvious trends such as those resulting from secular or 
seasonal variations could be subtracted from the series, the 
remaining series would contain random and independent data, 
and aeneras some of the not so obvious trends or patterned 
variations. | 


Many of the works in separating these "hidden periodicities 


were instituted by La Grange,?? Sheresaies, Whittaker,’ 


pokes, * 


and Schustenmee In all of these works, the random . 
independent series was assumed to be unimportant, and attempts’ 
were made to mathematically generate similar functions as 
those creating une Nnidden rrendar, | 

‘Schuster initiated the use of the periodogram, a method 
of pinpointing lesser periodic variations ina series. This 
approach was refined by such investigators as Whittaker and 
Pepinsonk” who showed that apparent major trends may only be 


a combination of multiple lesser trends av their Cormeen: 


peaks. 


46 4? 


It followed that Kolmogoroff , ‘- Weiner, and Cramér 


Should refine these attempts at separating periodic tendencies 





15 
by using Fourier series analysis to segregate obvious sine 
or cosine trends. All of these earlier theoretical works 
led to practical applications, particularly in the Teldec 
communications by ewes and his assistants at the Bell 
Telephone Laboratories. Applications in other fields have 
followed, until today, time series analysis and conmuniecenlen 
theory are major tools of the business and Be ononte world 


as well as the world of science and engineering. 








CHAPTER Sint 
HARMONIC ANALYSIS OF RANDOM FUNCTIONS 


A simple means of analyzing a function is through a 
resolution into its harmonics by the use of a Fourier 
series. Where the function has one independent variable 
and one dependent variable, f(x) = z, the one dimensional 
or single Fourier series iS appropriate. A two dimen- 
Sional function with two independent variables and the one 
dependent variable, f(x,y) =z, requires the use of the 
gdeuble Fourier Senne A discussion of these mathematical 
tools te contained in Appendixes A and B. Where random 
processes or random functions are being considered, the 
classical harmonic analysis is not directly applicable. 
However, it has been shown that a “type” of Fourier analyeue 
is still possible. This theory, which has .grown out of 
statistical theory of communication, provides the erie eee 
with additional tools to assist in analyzing random 
phenomena. 

Two very important properties of a signal or function, 
whitch for our purposes serve to characterize the function, 
are the power spectrum and the autocorrelation function. 
Prior to presenting the mathenatical background for cal- 
culating these properties, it may be well to offer a lay- 


men's definition of each. 
16 





EYE 

A. Power See Lat 

This is a term used in basic communication theory which 
has been borrowed from the field of elecUricivy aa aries 
analogies have been drawn to describe the difficult concept 
of spectral analysis, one of the better ones being that 
discussed in Granger and fates one The author will 
attempt to draw an analogy of similar nature. If a person 
stands on the corner of 42nd Street and Broadway, during a 
very busy semueae his ears will pick up many and varied 
sounds. Each of these sounds will appear to be different, 
although in many cases the "sound" which one hears is in 
meality a combination of many sounds coming from a range of ' 


_ ot 


Signal creating devices. 


BEES HATTER 


Figure,3.1 The Human Ear Performing Simple 
Spectal Analysis 
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The different sounds are a result of many signal 
generators, emitting different frequencies. The human 
brain is equipped with a device, not unlike a filter used 
in electrical circuits, which isolates each signal and 
attempts to analyze it. Sometimes it is successful and 
reports, "That is an automobile horn, and I must get out of 
the street," and sometimes the cacophony is so great that 
the individual signals can not be separated and interpreted. 
The brain is performing a type of spectral analysis when it 
selects, from a number of frequencies that the ear is hear- 
ing, individual signals and classifies them according to 
amplitude and frequency. Of course, in communication theory, 
and perhaps even in our analogy, there may be so many fre- 
quencies that it is only within our capability to identify 
and classify within ranges of frequencies. TLS ee ee 
often the case where a receiver is obtaining both intended 
or predetermined signals and "noise" which is added by some 
unknown and random cenenatons The mathematical process of 
decomposing a composite Signal into its component frequencies 
and determining the proportion of the total signal at each 
frequency has been called Power Spectrum Analysis, or it 
is often abbreviated as Spectral ered 

Spectral analysis has been used more and more in recent . 
years to analyze economic time series for the pacineee 
world. S\donmne - describes the use of this theory to separate 


the many factors which influence an economic series, such as 
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Br aaiiate and secular trends. Here, as in electrical and 
‘communication work, we discover the presence of “nome 
or factors which we cannot interpret as to origin, or 
predict as to future behavior. 
bee HUCOCOrRe abion Function 

The autocorrelation function may be described as the 
new function created when a function is correlated with 
itself. However, in order to provide Q more meaningful and 
definitive picture, we must look at the definition shown in 
Ekeblaa. 2 This text points out that the degree of auto- 
correlation in'’a time series indicates how well we may 
predict some as yet unknown future point in the series from 
the data available. A very Simple analogy can be given. If 
the sales in a certain business fluctuate from month to 
month as shown in Figure 3.2. and we compare each month's 
Sales with the preceding month's sales, we ae see that fol-_ 


lowing every high sales month there is a lower sales month, 


and vice versa. 


er 
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Figure 3.2 Sales Function 
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It is very simple to predict that in the future for every 
high month, there will be a lower month to follow. However, 
if we increase our "lag" and compare only every second month, 
we find that it is harder to make such predictions. If we 
go further and compare only the third months, we lose even 
more Penieidence in our predictions. This is a very 
Simplified explanation of Figure 3.3 whitch generally depicts 


the autocorrelation function. 


© \ 2 3 4 Ss & 4 ce => 
Lag in months 


Figure 3.3 Autocorrelation Function 


Figure 3.3 shows that: the autocorrelation function is 
maximum at a minimum lag and generally decreases as the lag 
is increased. Relative maxima in the function appear which 
Serrespond to basic cycles existing in the data, but unless 
the function is periodic, i.e., repeats itself identically 
at some latter period, the value of the autocorrelation 
function will never approach the value determined at Zero 
lag. Likewise, as the lag increases, our confidence of 


prediction decreases. 
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While we will not make use of spectral analysis or the 
autocorrelation. theory to separate and analyze signals, make 
predictions, or even use the correlation properties of the 
latter, we will attempt to take advantage of the properties 
of each to characterize or represent functions which are 
generated by our porous-media samples. For axaniouen the 
autocorrelation function has a property of retaining all of 
the harmonics of a given function, containing no new ones, 
and discarding the phase angle. The usefulness of this 
Peoperty shall become apparent as we explain tne applicavion 
of these statistical methods to the work of characterizing 
porous media. 

With the background definitions of the power spectrum 
and the autocorrelation theory, we may now show the mathe- 


matical derivation of each. 


weer eriodic Functions 
Lee?“ shows that a periodic function may be expressed 


as a Fourier series, as discussed in Appendix A. 





en 
ix ee AG. oe a cos 22% +p sin SY (aan) 
— > n L n L 
2 6 iar ; 
where: 
Tye 
aoe nix 
2 7 fi f(x) cos 7 dx a2) 


~L/2 





Ze 


yy Ae 
i fee aa ax | sua) 
-L/2 


oO 
i 
md 10°) 


L = Fundamental period of the function. 
x = value of the independent variable at equally — 
| spaced points along the function, such that x goes 
eigen, ly (We) Mes 
n = number of harmonic terms to which the series is 
arenncd | 


Equation (3.1) may be rewritten: 


n =o ; X 
ti.) = F(n)e J"T (3.4) 
| woe | 
im whiich: 
Eine ere ale)e J for ees; 1, Zoe ere 
Za 


j “= the imaginary operator, equal to the square root 
of minus (negative) one. (Some texts use the 
notation "i," and a more definite description may 
be found in any text dealing with complex variables.) 


Combining (3.2) and | Jj.) 
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Equation (3.6) is the Fourier Transform of the periodic 
function. Lee shows that this equation represents an 

"analysis" of the Aoi tudes and phase angles of each sinusoid 
into which function, f(x), is resolved. F(n) is also called 


the complex spectrum for this reason. 
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To separate the amplitude and phase characteristics of 


the complex spectrum, F(n), we may write it as: 
. eae | ey Aas 
F(n) 5 NAL vi exXp j tan Cz) SiR.) 
n 


where the amplitude is depicted by | F(n)| : 


DO 


[F(n)| = 1 a+b, © (3.8) 


and the term 8, indicates the phase angle: 


nN 
= a6 % | | 
of me (- 3s ) (3.9) 
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Combining the above, Equation (3.4) may be written: 
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2 a, *d, Cos a °, Grae) 


To analyze possible correlation between two functions, 
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we use the function, ne 1.2° where: 
1 +L/2 
— es + 2 
te 1,2 (ye) de ee f, (x) f,(x*T) dx (oan 


We call the left-hand member Ol Baa tonms an the 
correlation function at a given lag (7T). By multiplying one 
function by the other, at various values of VY , we can deter- 
mine the points of best correlation. 

If we adopt the special condition, where f,(x) et tx), 


we can write (3.11) as: 


ene 
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Lee shows that the autocorrelation function is also 


Sauay bor 


= . x 
. lPn)| 7 e nny T (3.13) 
n==-0 


When YT equals zero, Equation (3.13) reduces to: 








Pie : . ©o 
A 1,1 (“ane < 7 s ee (xdxe v PF, (n) : (3.14) 
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or the mean square value of the function equals the sum, for 
all harmonics, of the square of the absolute value of the 


spectrum. The power spectrum is then defined by: 


b (n) + JF, (n)| * | - (3.15). 


It can be demonstrated that the power spectrum and the 
autocorrelation function are Fourier transforms of each other. 
As previously stated, all functions consisting of identical 
SinuSoidal components have the same autocorrelation function, 
regardless of phase angle. Thus, the power Spectrum, as a 
Fourier transform of the autocorrelation function, is like- 
Wise independent of initial phase angle and having determined 
the latter, eye former is uniquely determined. 


Lee also shows that these functions may be written in 


terms of Fourier coefficients: 
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then: 
é 6 > 2 
$ (n) = [F(n)]* = 3n *Pn (3.16) 
and: | 
Rc ge Zee nn 
Say (TT ea bs *by ) cos Ser ete) 


From Bquations (3.16) and (3.17), it is@pessiolemuc 
calculate these desired properties of a function with the use 


of a computer program. (See Appendix C.) 


Wee Aperiodic Functions 

Whereas the tool of analysis for handling a periodic 
function is the Fourier expansion, it is necessary to rely 
upon the Fourier integral to analyze the aperiodic function, 
The derivation of this integral and the steps necessary to 
evaluate the power spectrum and autocorrelation function for 
this Brae of function are not required within the scope of 
this paper. A discussion of the theory is contained im) texgse 


D3 pal 


macn as Tolstov or Lee. 


EB. Random Functions ° 
We define a random function as follows: 
1. The generating process cannot be completely evaluated 
and a precise and unique determination of the value 
of the dependent variable for any value of the 


independent variable cannot be made. 
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2. As a consequence of (1), the value of the dependent 

variable cannot be predicted at any given future 
pPOLnv= | neveme. 

In other words, a random function is one that is generated 
by some unknown source which dictates the manner in which the 
function behaves, but, due to its stochastic or probabilistic 
nature, cannot provide the actual value of the function at 
any given moment in Ble: An example of a random function 
might be the noise that a radio picks up in addition to the 
intended message or program. ‘This noise may be the results 
of thermal emissions, statie, or spillover from other fre- 
quencies, but an exact description of the noise is not pos- 
Sible because of the complexity and random nature of the 
merse source. It might be noted that in most economic 
series, although containing obvious trends caused by seasonal 
em secular variations which may be isolated, Certain strong 
random GOnponenvs exisb. 

Because of the problems created by this "noise," it is 
of great import in communications work that analysis of such 
functions be a possibility. In order to analyze random 
functions, we find it necessary to introduce certain re- 
strictions, and assumptions. First, we assume that the 
function has.an infinite period and that during this period, 
the function will trace a pattern, such that the value of 
the dependent variable will pass through, or infinitely | 


close to, all possible values which are consistent with 
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the “ergotic hypothesis" as espoused by Maxwell, and 
discussed in Lee. Accepting this theorem, we may consider 
our random generating process as generating a group of 
finite length random functions rather than one infinitely 
long function having the combined characteristics of the 
entire grouping. The individual functions follow the nature 
of the infinite function on an average value basis and ina 
probabilistic manner. By the latter, it is meant that such 
properties as mean value crossings, or Zero value crossings 
per unit length Follow the probability distribution of the 
parent function. These finite functions will be referred to 
as an ensemble of random functions. They are all generated 
by the same source. This grouping of the ensemble is a very 
necessary step in the analysis of random functions, as 
realized and set forth by Gibbs.°? 

In such an ensemble of functions, each of the finite 
functions will be called member functions. It is important 
to realize that we assume that each of these member functions 
is generated by the same source, although we cannot mathe- 
matically formulate the operation of the source. This being 
SO, we sense that properties which characterize any member 
of the ensemble, Lyd haan characterize all re as well as 
eae infinite function from which the ensemble was created. 
One of these properties would be the autocorrelation function, 


ba 
€4, (Tee M™ Ls ti(ee, (am ax (3.18) 
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which we assume exists for all values of Y. Another such 


56 developed the 


property is the power spectrum. Wiener 
theorem that the autocorrelation function of a random func- 
tion and the power density spectrum are related to each 
other as Fourier Cosine transforms. The fundamental assum- 
ption underlying this theorem is that the autocorrelation 
function exists and is common for all member functions of 
the ensemble. Bie nmeneena Tukey?” State that, although 
Equation (3.18) is frequently called the autocorrelation 
punction, a more appropriate name for the term defined by 
Equation (3.18) is the autocovariance function. They state 
that the autocorrelation function should be defined as the 
normalized autocovariance function obtained as the ratio 

of the latter auncenen at any lag to the aurerion AGeacig@ 
lag. This affords the autocorrelation function the property 
of having the value of unity at zero lag. However, it is 
convenient for our purposes, and entirely consistent with 
the definitions, to use the names interchangeably for 
characterizing porous media. 

In order that convenient equations might be developed 
which will allow us to calculate these properties with the 
aid of a Aone, we refer to Granger and Hatanaka??: They 
Show that for.real processes the autocovariance,4,, may be 
written: 


Map : Ae cos Tw f(w)dw : (3.19) 


Ti 
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Formal inversion provides the power spectrum 


oe = Fn (4, + oat cos sw) OB 200 


J 
For finite data, using numerical integration, an 


approximation to Equation (3.19) is: 


° 


. n=T a Ks 
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where n * number of data points 


The estimate of Equation (3.20) becomes 
aw 1 n-l 
ae) ew omen uence Os i ) | (3.22) 


memnane” shows that this is a poor estimate, but that 
it may be improved by properly weighting the values of the 
autocovariance terms, so that 


aA 1 m--] 
art a = oT (eorota) 5 es Foot ee) cos J) (3423) 


where "m" is arbitrarily selected by the user and de- 
Moves the number -of lags used. It should be less than 


a 


one-third the value of "n." h,(w) represents the value of 


the weighting factors. 
There is considerable discussion concerning the Choices 
of the weighting factors, A(w). However, the two most 


widely used, developed by J. W. Hamming and J. Von Hann, 
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are very similar and for this paper we will use the latter's 
system, now commonly known as the Tukey-Hanning estimate. 


It has the values 


h #* 5 E + cos nn (3.24) 


The actual equations for estimating the covariance and 


the power spectrum may now be written 
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The Es terms are the raw estimate of the spectrum, and 
the ue terms are the smoothed values. The smoothing 
operation of Equation (3.27) is apparent from the form of 
the equation. ° 

To illustrate Equations (3.25) and (3.26), we might 


examine two very simple functions. Function A is an on-off 


function as shown in Figure 3.4. 





uk 





Figure 3.4 Example Function A 


By examining the figure, it is seen that the function 
repeats every. eight units and that the autocovariance 


function has the shape shown in Figure 3.5. 
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Figure 3.5 Autocovariance Estimate of Function A 


Referring to Equation (3.26), it can be seen that the 
power spectrum will have a peak value when the frequency 


parameter "j" is such that the cosine term, cos flee , has a 


period coincident with the period of the autocovariance 
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function. As an example, taking m = 150 lags, the power 
Spectrum reaches its peak value when the frequency parameter 


i is equal to 37.5. The power spectrum appears as shown 


in Figure 3.6. 





0 = o_o eee” 
° 10 2.0 30 


_ 


Figure 3.6 Power Spectrum of Function A at 150 Lags 


These results are easily checked by observing that the 
original function passes through 37m cycles in 150 units. 

A second function which can be used to demonstrate the 
nature of the analysis performed by use of the auto- 
correlation function and the power spectrum is Function B. 
Function B is a combination of two simple on-off functions 
which have been superimposed. The frequency of the first 
function is four cycles per eighty points. The frequency of 
the second function is ten cycles per eighty points. 
Function B, as shown in Figure 3.7, has very little apparent 


regularity. 
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Figure 3.7 Example Function B 


The power spectrum for 50 lags shows peaks at 5 and 
12.5, which were the individual peaks for each basic fre- 
quency used in the combined function. In other words, 
Equation (3.26) effectively selects the fundamental fre- 
quencies which constitute a function, even though they may 
not be apparent individually. Further, a measure of the 
relative contribution of each frequency is offered. 

The computer program "COVAR," discussed in Appendix C, 
has been written specifically to perform the mathematical 
operations set forth in Equations (3.25) and (3.26). The 
above figures were made from data and results of that pro- 


gram, run on Functions A and B. 








CHAPTER IV 
EXPERIMENTAL PROCEDURE 


A. General 

In this thesis an attempt is being made to classify 
or characterize porous media. It is our goal to differ- 
entiate between media which are not statistically identical 
through the siete of vatfious forms of statistical analysis. 
In order ‘to apply the various tools outlined in Chapter III, 
we require a variety of samples which we know to be Sictimi lar. 
In this chapter, we will discuss the selection and preparation 


of these samples and the collection of data therefrom. 


8 


B. Preparation of Senet ae 

Data used.in this thesis were taken from two types of 
porous-media samples. The first sample was a photomicrograph 
of an actual porous medium. The other samples used were 
mee oared by the author by imbedding multi-size spherical 
glass beads in Deis Wer 

The photomicrograph is a photograph of an actual porous 
medium taken through a microscope. In the sample used, the 
photographic image wa8 enlarged ninety-two times from an 
original size of approximately .12 centimeters by .16 
centimeters. This photograph was taken from the paper by 

60 


Fara and Scheidegger. The original medium was sandstone 


With a porosity of 26.3 percent. . 
34 
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The other samples used were created by the author in the 
following manner: 
Materials: 
(1) Laboratory glass beads, three sizes: 
ae 4mm, Smm. 
(2) Natco Premium Quality Clear Polyester Casting 
Resin, NL-600. 
(3) Violet dye (water base). 
(4) Waxed paper cups. 
Procedure: 
7 (IB) lakes dye and resin to achieve a colored, 
Slightly opaque plastic. 
(2) Place liquid plastic in waxed paper container. 
(3) Drop glass beads, in varying proportions as 
required, into the plastic. 
(4) Place ‘in over for one day (30 degrees centi- 
grade). Allow to set. 
(5) Cut through sample in a random manner to 
realize an unbiased cross sectional surface 
of the mold. This surface area represents a 
model of a porous medium. Where the glass 
beads are apparent is taken as solid-space, 


and where plastic is apparent is taken to be 


pore space. 
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For this paper, three samples were prepared. Sample X 
consisted of 3mm. beads only. Sample Y was a mixture of 
equal volumes of 4mm. and 5mm. beads. Sample Z was a mixture 
of equal volumes of 3mm., 4mm., and 5mm. beads. The cut | 


Samples were: approximately two inches in diameter. 


C. Preparation of Data 

Following the preparation of the samples, photographs 
of each were made by the Photo and Graphic Arts Department, 
University of Kansas (see Appendix D). The images of the 
created samples were expanded six times, while that of the 
Scheidegger photomicrograph was already of appropriate 
dimensions. An 80 x 80 grid, with ten units per inch was 
overlaid on each photograph. At each grid point, the data 
in point was recorded as being a *1 if a pore space (plastic), 
and a -l1 if a,solid space (beads). These 80 x 80 data 
points were easily converted to input information on eighty 
aogier computer cards, for use with the programs written 1 Oi 
this thesis. (The negative values were read into the 
computer as fixed point zeros and changed internally to their 


correct values as negative ones.) 





CHAPTER V 
DISCUSSION OF RESULTS 


The mission of this project was to provide means of 
statistically analyzing porous-media, to show statistical 
variation in media known to be different, and to offer 
direction and inspiration to further work in this area. The 
author fully realizes that from the limited number of sam- 
ples dealt with in this thesis, there is insufficient evidence 
to prove that the far reaching goal of statistically charac- 
terizing all media can be accomplished. However, in this 
chapter, we will offer indications that the above outlined 
missions have been successfully completed and that success 
of the overall goal mentioned is entirely possible. It is 
hoped that from the results shown, others will be stimulated 
to continue research toward that end. 

Two basic computational tools were used to statisti- 
cally analyze the data taken from the porous-media samples. 
Although discussed in Chapter IV, it might be well to point 
out again that these data were taken at equally spaced points 
along lines randomly drawn across the cut-face of the samples. 
Where the line passed through pore space, a value of *1 was 
assigned to the function, and where the line bisected solid, 
@ -l1 value was assigned. A series of these points taken in 


W 





38 
order along the line then formed a simple orthogonal function 
with the +1 to -l cycle arranged in a seemingly random pat- 
tern. Based upon the as sunpeon that porous media are 
Statistically constituted, it was hoped that several of these 
lines or functions taken from the same medium would be 
statistically similar while those taken from different media 
would vary by what is known as a "significant difference." 

The first method used to statistically analyZe the 
functions was to fit each with a single Fourier series using 
"1D Fourier" computer program. (See Appendixes A and C.) 
Although each function was considered as being random in 
nature, it was analyzed over a finite length. The Fourier 
coefficients were determined and then combined according to 
the Lee estimate of the power spectrum, one > a* ¥ le SO 
that the initial starting point or phase angle might be dis- 
regarded. These values of Ci» as eee of the power 
spectrum, are measurements of the contributions of the 
fundamentals frequencies constituting each function. ‘This 
is essentially the technique used in the paper by Fara and 
foneneeerer.°! | : 

Typical results of the single Fourier series approach 
are displayed in Table 5.1 Shown are the CL values for four 
lines across the same media. As can be readily observed, 
there is very little consistency among the values, indi- 
cating that this method yields a poor estimate of the power 


oo. 


spectrum. As a result, the simple Fourier approach was not 
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pursued to any great detail, although later in this chapter 
we will compare the results of using this approach with those 
of our second technique on an identical sample. 

Because of the poor results achieved in the single 
Fourier approach, no analysis of the two dimension cross 
sectional faces of the media was attempted. However, a 
computer program "2D FOURIER" to determine double Fourier 
series coefficients for a surface was developed, tested, and 
is fully operational. (See Appendixes B and C.) 

The second statistical tool used to analyze the data 
was the determination of the communication theory properties 
mR ie functions or lines. These properties, the autoco- 
Variance function and the power spectrum, are discussed ia 
Chapter III. The properties were calculated using Fortran 
equivalents of Equations (3.25) and (3.26) in the computer 
program "COVAR." (See Appendix C.) As a part of that pro- 
gram, the power spectrum was plotted for each of several 
lines per media type, as well as the 90 and 95 percent 
confidence bands. These bands were created by multiplying 
the smoothed power spectrum by constants presented in Granger 
and Hatanaka.°@ The constant used is a function of the 
"equivalent degrees of freedom," an approximation chosen to 
allow for the fact that. the series under consideration is 
not necessarily exactly normally distributed. These bands 
indicate 90 or 95 percent confidence that the true value of 


the power spectrum lies within that range of values. In 
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other words, 90 or 95 percent, as the case may be, of all 
power spectra values from the Same population should lie 
within that range. Generally, the approach was to compare 
data lines from the same media to see if the observed plots 
of the power spectra were similar. Then it was determined 
whether the values of each function within the sample, at 
each frequency, fell within the confidence bands of the 
other lines within that sample or type of media. A similar 
technique was used to compare the power spectra profile of 
one sample media with that of another media. 

For these comparisons, the effect of sample lengths was 
investigated. oSamples of various lengths were created by 
combining the eighty point crossings developed in accordance 
with Chapter IV. The crossings were shuffled to reduce 
bias and eliminate peaks which would occur in the autoco- 
variance function at eighty point intervals. These peaks 
would be due to the similarity between adjacent lines drawn 
across a porous media. It was discovered after a series of 
trials that the longer the lines, the less bias there was 
to the idiosyncrasies Of (SINS erweroOssi NaS. 1.6. mine 
might pass through continuous solid or pore space for an 
many as twenty or thirty points without interruption or 
@maange. Longer lines tend to minimize the effect of these 
rare exceptions. As a result, lines of 3200 data points or 
40 crossings were utilized, this being the maximum number 
of points that could be handled on the computer at one time. 


It might be noted that lines of 1600 data points afforded 
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very similar results, indicating that the lines of 3200 points 
Pere Of SULTAIC1eny length. 

It was discovered that the number of lags selected, or 
more properly stated, the length of the maximum lag selected, 
played a very important part in the aneieicee Again after 
many trials, it became apparent that when the maximum lag 
selected was too short, the apparent frequencies were 
"packed" at the very low power spectrum values. In other 
words, with many frequencies constituting the function, 
eocn band or erouping of frequencies was too wide or contained : 
too many frequencies to effectively delineate the function. 

To long a maximum lag had the erfect of causing random spikes 
in the power spectrum which were characteristic of an in- 
dividual line‘but not of the media as a whole. Theses *epinera 
might EE a result of sampling variability. Also, 
at long lags, values of the autocovariance function vended 
to distort the power spectrum profile. In light of the. 
above, a lag of 40 was selected as Senne optimum. It had 
the effect of smoothing the data and still offering good 

Me orimination, | 

Figures 5.1 ene 5.3 offer three plots of the power 
spectra of different media. Appendix E agreeing the com- 
puver results from which these lines were plovledsandmiume a 
computer plots. The three profiles are obviously different, 


but it is important to discuss these differences in a more 


absolute statistical manner. As Shown in the plots, the 





Ne 

95 percent confidence band of one of the lines is Super- 
imposed upon the PM yaNies By looking at each figure, it can 
be seen that generally the power spectrum of as OUNGeaor 
Pouter" bee lies within that confidence band. More 
specifically, the outer line is within the band at 24 points 
out of 40 in Figure 5.1. This is the equivalent of 60 per 
cent of the time. In Figure 5.2 the result is 23 out of 40 
(57.5%) and in 5.3 it is again 23 out of 4O (57.5%). 

While these results violate the 95 percent criteria, 
it is important to note that many of the points are very 
evese to the limits Sra as only two examples were displayed 
per sample, variability of samples has a great effect. 
idleeaaley many more lines per sample should be obtained, 
pleuted and analyzed. Then the confidence bands around whe 
mean of those lines should be superimposed upon a plot of 
the lines. Looking at all of the power spectra values that 
lie within these bands, the 95 per cent criteria would be 
epproached. Using this method on four lines of shorter 
length from these same media, results very close to that 
Griveria were realized and as more samples are considered, 
it is expected that the 95 percent criteria would be ful- 
filled. However, this was not proven in this work. By the 
time the tools to analyze media had been developed, and made 
available in the form of fully tested computer applications, 
and when the results discussed above became apparent, there 


was insufficient time to create or acquire additional samples 


for this thesis. 
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Cross comparing samples, there is very little correla- 
tion in comparison to the correlation within the same media. 
Between Figures 5.1 and 5.3, which are very similar media 
(see Chapter IV), the base lines are within the other's 
confidence bands only 12 of 40 points, or 30 percent of 
the times. Between Figures 5.1 and 5.2, the result is only 
one out of 40 (2.5%), and between 5.2 and 5.3, it is one out 
of 40 (2.5%). Visually comparing the different spectra, it 
can be seen that the 3mm. beads sample, which only has small 
beads in its make-up, has lower values than the other 
samples at low frequencies, and higher values than the others 
at higher frequencies. This is to be expected as there is 
a greater fluctuation in the +1 to -1 cycle in the original 
3 mm. bead data. There is less opportunity for long strings 
of +1 or -1 values due to this composition of the sample. 
Similar analysis can be made between the other samples. 

It is felt then, that this is a definite indication that 
the general shape of the power spectrum shown is indicative 
of the composition of the media from which the data is 
taken. Spectra from different media are not only visibly 
different, but also it is noted that the differences appear 
to be statistically significant. Again it is repeated that 
many more samples of various types need to be analyZed to 
prove this conclusively, but the technique is readily avail- 
able, and indications: are that the results would be 


favorable. 
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Of interest, Table 5.2 and Figure 5.4 display the com- 
puted autocovariance estimate and the power spectrum of the 
actual rock sample discussed in Chapter IV. 

A brief comparison was made between the two techniques 
used to analyze the functions, the single Fourier series and 
the power spectrum. Figure 5.5 shows this comparison. The 
function analyzed was the line drawn across the porous media 
shown in the Fara-Scheidegger paper. (See Chapter IV). As 
can readily be observed, there is very little consistency 
in the results. This is further indication enae the single 
Fourier series approach is not particularly suited for this 


type of analysis. 
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POWER SPECTRUM ESTIMATION 
FOR DISCRETE TIME SERIES 
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C(O) is an indication of the porosity of the media. The 
actual porosity of the sample was 35.65 percent. 
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CHAPTER VI 
CONCLUSIONS AND RECOMMENDATIONS 


In this thesis,\an attempt has been made to provide 
effective methods of State vice ana aZane: porous media. 

It has been shown that the properties derived from communi 
Mion theory, the autocovariance estimate and the power 
spectrum, are capable of differentiating between different 
media. Further, it would appear from the work with the 
simple representations of media which were created for this 
problem, that real porous media lends itself to this type 
©f analysis. 

Perhaps the most significant conclusion that can be 
drawn from this thesis then, is that the techniques of 
statistical communication theory are adaptable for this form 
of analysis. Although the classical Fourier approach did 
not provide consistent results, the Tukey Hanning 63 power 
spectrum estimate offered good discrimination between the | 
Gifferent media. 

There are several obvious steps required to further the 
proof and development of the theory presented herein. A 
large number of samples from each simple type of media should 
be analyzed to affirm the consistency of results. Then new 


samples should be fabricated with varying proportions of the 


De 





Be, 
different size beads, to determine how fine a discrimination 
can be made. | | 

If these tools continue to offer good results with 
idealized media, the next step would be to test them on real 
meck Samples. ~such Statistical methods as Discrimination 
Function Analysis or Factor Analysis might be applied to the 
data to improve the differentiation between "different" 
Medic. | 

Finally, if these techniques prove successful upon real 
media, it would be logical to then attempt to relate the 
statistical properties, such as the power spectrum, to the 


4 


static and dynamic fluid behavior of the media. 
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APPENDIX A 
SINGLE FOURIER SERIES 


In many technical problems, iteis necessary bo deal wing 
complicated periodic functions, that is functions Pe 
F(x + T) for all values of x, T being the period of the func- 
tion. Use of the Fourier series allows us to write these 
complex functions in terms of simple trignometric functions. 
Although Pourwer series applications are not. restrre ved sre 
periodic functions, as explained in a Speci a ieeame 
must be taken in applying the techniques, the discussion of 
which is beyond the scope of this paper. Most texts on the 
Wse Of the Fourier series provide the necessary instruction 
mietuhis area, €.£., Move rou 

Although the complete theory of Fourier series expansions 
are covered in detail in such reference material as Tolstov, 
Gaske11, 6° Hildebrand, ©7 and others, a resume of the simple 
derivation for a periodic orthogonal function is ideacem 
herein to assist the reader in the basic understanding of the 
Mses of the harmonic analysis that follow. ~ 

If we have a function, such that there is one independent 
variable and one dependent variable, taken between May Wuntine 
=i = 0 < L, we could write the funetren: 

- +7, 


Z, E(x) or; Ze S xdx 
-L 
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oz 


We might illustrate the same function: 





The data set, (x, , 


Z4) is such that x eee is (X44 — Xe 
constant and x, ~ *L, the value of "k" being even. 

We then assume that an equation using simple sine and 
cosine terms may be written, which will represent the func- 


Pron. as follows: 


& et ee TX 4 A Oe nx 
vee f(x) - a, cos “7; + a, cos 4% . . . &, cos “FP 
+ yee in SRK 4 ny DRX 
by sin be, sin => eo db sin “> 
onae 
Z* f(x) «2+ 5 (a. cos M*® +p» sin BY) (1) 
a, n L n L 


The degree to which this equation, so written, will converge 
upon the function, may be shown to be dependent upon two 
maecors. The first factor is the Aitber of harmonics (n) to 
which the series is expanded, and the second is the number of 
data points observed in relented to the basic length. The 


author has achieved an accurate fit to the degree of .974 on 


en = 
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an eighty point function of length 2Mu, with data values 
manging from =I Go 71) The series was expanded emer... 
harmonics. By holding either of the stated factors constant, 
and decreasing the etner, a Voss i accuracy of represent. suem 
was realized. The reason for the "2" in the denominator of 
the first term of equation (1) will become apparent as we 
find the coefficients 2, and I. of the series. By multiplying 
each side of the above equation by cos nx, and integrating 
between the limits of *L and -L, Gaskell shows that the 


equation may be rewritten: 


et, ag 


nix - 2 nx a 
. f(x) ces = ax a. ia cos” “T= ax La, 
or: 
+1, : 
ee Scuola. = dx (2) 
L -L . 


With the a, term as the first term of the equation (iB) 


utilizing the same steps, we discover: 


wt al +7, 
J f( <jmcos Lo}mx a, = f f(x) dx = f tin dx 2 La, 
-L -[, -{ 
aa call 
a cis iQ ed x Sa 4 S ii ees (o)mx ax (oy) 
O L ot ls 1 L 


a 
Gaskell further notes that oa is the average value of 


f(x) over the entire ine 2 eens the trapezoidal rule to 


numerically integrate Equation (2) over its limits, and 
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js 


utilizing equally spaced points (Ax = ) whence inne Tamar 


= 


. . . ky, we can approximate Equation (2) as follows: 





; + 1 *=*k-L 
a. - fa » AX 26 a onus ne? Z cos Di - 
2 ll L 
Woe 
N26. Boe cos’ a ec ocn (4) 
k 2 ile L 


It is shown by Whittaker and Ropinaenne, that use of the 
trapezoidal rule produces coefficients which satisfy the 
least squares criterion of a "best" fit to the existing data, 


mat 1S, 
k-1 


i a | 


| (Z, = Px, )° = Minimum (5) 
re) 
where PX, is the polynomial approximation to be observed Zee 
The same reasoning as above can be enlisted to determine 
fee coetficient by» except that the sine terms ane AEN AG) ENG 


the end points, (sine (0) = O) and b, becomes: 


1=k-1 


oy z. sin OT%; (6) 





Now, assuming we have discrete data Z. for each value 
of Xs taken at equally spaced points between the limits, we 
can now compute the values of the coefficients, and having 
so done, substitute back into Equation (1), which becomes a 


one-equation representation of the function. The coefficients, 


a. and re serve to characterize, in one way, the function. 





65 


However, this means of characterization has the distinct 
disadvantage of being dependent upon the phase angie of the 
Sinusoids into which the function is resolved. 9 AN aints 
disadvantage may be overcome in several ways. Some of one 
more complex characterization techniques, which we shall 
discuss and compare in latter sections, are the use of the 
autocorrelation function, and the power spectrum. These 
functions are non-dependent of phase angle, as we oil 
later show. A less complex technique, which still dismisses 
the dependency upon the nasceenee is discussed in a paper 
on wind generated sea waves, distributed by New York 
University.?2 An estimate of the power spectrum is made 
in the following manner, using the Fourier coefficients: 


Estimated Power Spectrum cf bs ae i db“ (7) 


Although this estimate is a characterization of the 
function, and does dismiss the dependency upon the phase 
angle, the same reference notes that the quantity ae is a 
very unstable estimate of the energy at any frequency. | Only 
through a weighting process of the Ch values, would an 
accurate spectrum, such as that determined by the Tukey 
method“ be recovered. Lee? uses a Similar estimate of the 


power spectrum, with the exception that he divides the above 


Sestimate by a facvor of Tour. 








26 
In that the basis for our discussion dealing with auto- 
correlation and power spectrum analysis evolves from this 
same reference, and since thes@ivision does NOL 4D ReCeiaaeme 
term for our purposes, we have used the "Lee" estimate in 
our work. This usage will allow a means of comparing this 
Simple estimate of the power spectrum with the values 


computed by more complex techniques. 





APPENDIX B 
DOUBLE FOURIER SERIES 


As in the case of a function with one independent 
variable discussed in an earlier section, it is also possi- 
ble to represent with a Fourier series, a dependent 
variable as the function of two independent variables. 

ve Ne) (1) 
A Poor ous discussion regarding orthogonal systems in two 
variables may be found in Tolstov.?4 

Thus, a parameter such as Zz may be considered to be 
oscillatory in two mutually pencenaicwlar directions, such 
as east-west and north-south, allowing Equation (1) to. 
lmepreesent the areal variation in z. If we regen: a 
function z to have a fundamental period of 2L along the x 
coordinate, and 2H along the y axis, then following the same 
reasoning discussed earlier in the Single dimension case, we 


Can write Double Fourier series function: 


m=M n=N | 
Z = Pp =! a) r rN a cos Timxk cos Tn 
XsY =o nzo)0 OM | On La: ais : 
oe o Sin tmx cos Tiny + C, e cos Tmx Sin Tiny 34 
’ L H g Js H 
d sin tmx sin tny (20) 
m,n a A 
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where 

IN = 4 ies=n =o 

Wal 1/ 

h <2 Ae ml = ©,0-e some of. © 

m,n 

r ia ee m>o,n>o 

met , 
M,N number of harmonics considered 


(Should be equal to or less than one 
half the number of data points con- 
Sidered in the x and y directions 
respectively, 79 
The reasoning behind the quarter and half terms follows also 
from the one dimension analysis, in that the end terms con- 
sidered in a numerical integration are:-averaged, according 
to the theory of the trapezoidal rule. 7© In the two dimen- 
Sion case, the averaging process must be carried in each 
direction, thus the quarter terms. Again, the problem is to 


evaluate the coefficients. Integrating Equation (2), the 


eeerricients are determinable from the following reiations: 
+H 41, | 
=i mmx Ty : 
Qnon ” TH - a i (xX, ye moos Tr COS —y* dydx (3) 
1 HO fb TI TT 
o mx ny 2. 
aes 7 - i i ( Xgy ese “pcos “ dydx (4) 
. 1 "HO SL TI 
C = i- f f f(x,y) cos ™ sin — ayax Cop 
m,n LH CL L H 
1 , H 7b TI TT 
« mx ny 
1 ra f f(x,y) sin == sin =)" dydx (6) 


-H ae 
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IG is possipie yuo numerically integrate Equations 3-6 
utilizing the trapezoidal rule. Again it must be remembered 
that the integration must be carried out in two dimensions, 
and special care must be taken with regard to the corner 
points. In order to clarify the following derivation, 


Figure 1 displays the nomenclature to be used. 
j-4 


y Um ees 950 Ba! 
j Z(i, 3) | 
18 © gle eres ene 


--J™=o 


gal fenvetsy IL 


According to the proper use of the trapezoidal rule, 
it can be shown that. for the one dimension case (one 


independent variable) 


where Ax = a 


Considering the second dimension of depth in our picture, 
when we are concerned with but one direction, we are able to 


write Equation (7) 
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+T, i=*k-1 


‘ a + 
i 7 (x ce | arom MM oan? a Y) Zz. (8) 
\ 2 1=1 15) 


However, in a specific case, where it is required that 


e 


integration be applied in two dimensions, as will be done for 


Equation (3), the following steps are taken: 


1 1 
a eG iy i Zs cloacley 
oH Zo SZ, ie Ke | 
og; Ax Ser SL = Zs 4 dy 
Loe = ee 2 = oJ 
1 gal 
a (A) 
, Axdy FOG : EKHon One Zia rt 
~ GLH 2 Z 
2 
‘ra \B)  fink-aggy | se) 
2 One) k,] +) 2% L@ ines + 
he) z 1 ee 
Sf i (9) 
j-l al 459 


where: (A) term that averages the corners in both 
Gdirectrons. 
(B) term that averages opposite vertical edges 
(less corners). 
(C) Term that averages opposite horizontal 
edges (less corners). 


(D) term that sums all interior points. 





71 


The number of points in either direction should be an 
odd number, requiring "“k" and " " to bé even numbers. Carry- 
‘ing the procedure to completion and reinserting the trigno- 
metric terms, which were Mire eu (6 simplify the explanation, 


we can write the coefficient as follows: 





+ + A 
a BS Lal 2b Bi “0,0 “kyo + fo... “kee _ “ka ve 
m,n _bit a. 2 2 cos ml cos nt 
| 2 
=Q_} (B) 
+; 2 GON i “k,j cos mm cos nny , 
i- jF1 o) a A 
g Riera 
+] 2: i om eae mT x . (C) 
i=] t D ee ©) ee cos nv 
ik $ 
jA- oe ie i (ep) 
z. . cos ™"*s Gog any 4 (10) 
=): te) ie Al . . 


Similarly, Equations (4-6) become: 


oo ee 
m,n 1k No (A) term in that on the corner 


inne Tl 
points, ic ll TODS) 5 


a 


b 


No B term due to the fact that 


xX and x, are equal to i and sine Bi 


term goes to zero. 


hae one 





ais COS 


a 
a 


ta = Zs 4 sin ™'%4 cos any «| | Cp) 
ee , H 


| js : pes bs 
non 
‘ad 
N 
O 
N 
a 
po 
n 
-~ 
aj 








V2 


No (A) term because sine terms go to] , 

zero. 

j= -1 (2) 
+ 

De Sone Zien § cos mm sin nny 5 + 

jF2 2 H 


No (C) term because sine terms go to} , 
Zero. 


jef-1 itk-1 eae ane (D) 
u a Z4 4 COs i sin j 
Sil sc Tie b - ie H 


No (A), (B), or (C) terms because sine 


Pemms FO tO Zero. : i 
| _ (D) 
As) oa Oe et oT) ; | 

s oF zy, sin mXs sin ONY 4 

Len si eil ’ ip H 





APPENDIX C 


COMPUTER PROGRAM TO DETERMINE COEFFICTENTS FOR SINGLE 
FOURIER SERIES (ONE DIMENSION) FIT OF ORTHOGONAL ‘FUNCTIONS 


Pur pe ce 


The purpose of this program is to determine the Fourier 
coefficients a, and b,, when a function, f(x), Ve Wal eeyel shy 
& Single or one-dimensional Fourier series. The secondary 
purpose is to calculate an estimate of the power spectrum 
of the function, which in turn is a characterization of the 
maaction. 


Language 


Fortran IV (IBM 7040 Computer) 


symbolic Dictionary 


Variable S/A* T/o** | Description 

x A I/O Independent variable data 
ipoin v.. 

Z A iy, © Dependent variable data 
pens. ; 

MZ A I Fixed point value for Z. 

A A O Fourier coefficient. 

B A O Poummer = CoOcii iC lem. 


= o--Single Variable; A--Array of Variables 
a re Variable; O--Output Variable 


(s 





Variable 


C 


CC 


BL 
KOW 
KOO 


Pir 


KKK 


ZO UM 
SUMCS 
SUMSIN 


SUMCOS 


AZERO 
I 
PMAX 


PMIN 


cae s 


Ann MH N 


WY Wt mu 


O © tonne 


Description 


Estimate of the function's 
power spectrum based upon the 
equation: 
Cc =A* 4B 
n n n 
Estimate of input function 
Z at each data point, created 
by using calculated 
coefficients. . 


2 


¢ 


@ 


Residual or difference between 
iipMietuUNnction Z and estimated 
function P, at each data point. 
A(N)© + B(N)* 

Number of data’ poimps: 
Half-length of sample. 
Number of samples. 


cample index. 


Mathematical expression input 
as being equal to 3.%4159265. 


Number of harmonics to which 
the Fourier series is to be 
expanded. 

sum of Z data point values. 


Sum of Sine and Cosine terms. 


Sum of Sine terms in equation 
to determine coefficients. 


sum of Cosine terms in equation 
to determine coefficients. 


First term of Fourier expansion. 
Data point index. 
Maximum value of variable P. 


Minimum value of variable P. 
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Variable S/A I/O , Description 
ZMAX Ss) O Maximum value of variable Z. 
ZMIN Ss) O Hannum value of variablews. 
RM AX S O Maximum value of variable R. 
RMIN S 0 Minimum value of variable B. 
N S ~- Harmonic index. 


Program Routine 

The program utilizes data points equally spaced along 
the function from a set origin, andebyauel he tt oOrleanmedua— 
valents of Equations: (4) and (6) from Appendix A, generates 
a Fourier series estimate of the function. The Fourier 
coefficients generated are presented as characterizing the 
function, and by combining them according to the “Lee” 
estimate, c(N)¢ = (A(N)* + B(N)©)/4, a Shateetenieine power 
spectrum which is independent of phase angle, is determined. 
When more than ‘one sample is run, the C(N) eee are grouped 
statistically for analysis. The program checks the values 
of the coefficients by creating an estimate of the original 
mane tion, using Equation (1) from Appendix A. A plot of the 
Peieinal function, Khe created estimate, and the residuals, 


or differences between the two is presented. 





MAAAAOAAY 


Po 


203 
204 
800 
801 


205 


2 


PROGRAM -~FOURIER- 


SINGLE FOURIER SERIES PROGRAM 


PROGRAMMED BY WD ALOENDERFER 


ONE DIMENSION FIT OF 3MM BEADS DATA 


400 POINTS 50 HARMONICS 


DIMENSION X(801)92(801) »MZ(801 


TAIT E GES ist.0 (30) sic | 
(0) 


LIDSH(32) » 
GALL CLOCK 
READ (5.102) (TITLE(T)sI=1913) 

FORMAT (13A6) - 

READ (591) Ks EL 

FORMAT (110910X9F10-0) 

WRITE(656000) 


ae en oe — me —2e 


- — oe we ee oe © - - 


AUGUST 1965 
(ONE DIMENSION) _ 


19A(30)98(30)0C(30930)sP(801)5 °° © 
30) »LR(30) sR(801) 


me me ee ee ee ee ee ee es ee ee ee we a ee ee ee 


— <= a ae ee ee ee ee ee ee oe ee ee ee i 


FORMAT (1H1s20Xe921HPOWER SPECTRUM VALUES) 


KOW = 16 
_ KOW IS THE NUMBER OF SAMPLES 
KOO i 
PIE an 
AK 
MK 
AMK 
KKK 


= 14159265 
— a y 
i ol | 
MK 
50 


KKK IS THE NUMBER OF HARMONICS _ 


DELX 
X(1) 
ZSUM 
SUMCS 0.0 
Z(1) Oe 
DO 203 IT 
me 1il) 
X(T) 
CONTINUE . 
READ(59204) (MZ(T)>5 
FORMAT (8011/8011/80 
DO 205 JJ = 1loK 
IF (MZ(JJ) - O 
Z( JJ) -1e0 

GO TO 205 

Z2( JJ) MZ (JJ) 
CONTINUE 

wor i2 IT 29MK 

Z2SUM ZSUM + Z(IT). 
CONT INUE 

SUME =(Z2(1) + Z(K))/26¢ 
AZERO 
DO 4 N 
SUMSIN 
SOMF 
SUMCOS 
FN N 
SUM= 020 


PIE/( (AK = %6")7/2ieq) 
~(PIE) 
020 


0 — 
I = 29K 
= 0.0 —_ 
X (I 


II-1) + DELX 


[=i ) 
lw 7S 
}8005800,801 


ee 


1 KKK 


QO ot 


0 
rae) 
= 020 


RN AO gD. - Sinn 


O1l/80l1/T1) 


1 — << os Om mtg ae te 


(2¢/AMK) *(ZSUM + SUME) 


me mt ee ee et ee ee ee ee ee ee ee ee ee —— — 


BEING TESTEDs 


.— a) oe ae a ems ets oie oe a or = a na ame ema = /—— a a oe —— = <= 


ete ee ed eh ee ee ee ee ei a i ee ee ee ee es ee ee ee et ce te ee ee ee ee ee ee ee et ee - = 

ee eee ee eee ee 

— ——-12 — Se he ce ce ee ee ee ee te ee eee ee ee ee ee ee ee — ow oe ee oe ee 

noe ae eS 
I ee 

— ie te i Sa EE ae et Oe Oe Oe 

ee re 


ee 


ee ee ee cme ce ee ce ee ee es ee ee ee ee es we we ee ee : : 


ALGEBRAIC INDEX EQUALS FORTRAN INDEX LESS ONEe 


DO 5 I 2 9MK 


SUMCOS = SUMCOS +  Z(1)*COS(I 


FN* PIE * X(1))/EL) 





NADNNAN 


687 


6005 


SUMSIN 
CONTINUE 


SOME 


_SUMSIN + Z(I)*SIN( (FP NE TE 3-960 1 Ae ee 


((Z(1) “$2 220) COS LENMR LE) 


—— =p =e ope oe ap = om Ow oe Ge = ae ee eo 2 ee ee ew Se ee eS ee SS Se ees ee ee ee ee See ee ew 2 ee Se ee See SP ea ap 2 ep ep 62 ep eS en S22 ae oe Se oe Ss op eee ee ee = a2 oe eo 


COMPUTING THE VALUE OF THE POWER SPECTRUM ESTIMATE e 
ESTIMATE 


C(KOOsN) 


A 


A(N) SQUARED + B(N) SQUARED DIVIDED BY FOUR _ 


2 Qed OD oe ep Ae GD Dee ek eS Gm een Ge GD ae eee eee + Gee Ae ee ee | 8 ee ee ee + 


2e/AMK * (SUMCOS + SOME) 
20/AMK * SUMSIN 
(N)* A(N) + B(N) 
CeVa0 


* B(N) 


CONT INUE 


WRITE(65687)(C(KOO9N) sN=1 950) 


FORMAT ((1Xs1l0F1004)) 
WRITE (656005) 


FORMAT(1HOsl1HNEXT SAMPLE) 


GO TO 686 


NMAX 
DO 8 


SUMCS 
DO 9 N=] KKK " 


FN 


SUMCS 


~ 1TB(N) 
9 CONTINUE 


DADA N 


214 


10 


tal) 


R(T) 


—— 


I 


N 
a 


Fs ae et ete = ead ee oe 


(K-1)/2 
1K . 
AZERO / 2e 





Se __ 


SUMGS HAIN) ee -GOS((EN* PIE © X(T) )/EU eee 
SINC(FN * PIE * X(1))/EL) 





“SUMCS 


P( 1) — mia) 


CONTINUE 
SO T.0.-2.1.3 


ae nik suki dandeaes teniaendaeaaaet 


me ee ee we ee rm ee em ce re ee ee ee en ee cr cm ee ce ee en em ee cr re re wm et ee ee te ee 


NE YE 


IBLANK 


IDASH 
my 10 | 
IDSH(I)= 


-13702990896 

27661634608 

-17997958192 

-818089008 
1932 

IDASH 


atl 
=> 
= 
= 


WRITE (65103) (TITLE( TIT) sIT=1913> 
103 FORMAT 


(1H113A6//917X»8HOBSERVED 932Xs10HCALCULATEDs 


132Xs8HRESIDUAL } 
WRITE 
1(IOSH( 1) sT=1932) 


SS er 


(69104) (IDSH(1)9I=1932)s(IDSH(1) s1l=1932)9 | 


Sa eS eet eh ee pen eee ee ee te ees te ee ee ee ee a es et es ee ee ee ge ee ee ee ee ee ey 


104 FORMAT (5X932A199X 9 32A1 99X 932A ) 


ZMAX = Z(1) _——— ll Gein ee... 
ZMIN = 2(1) 

PMAX = P(1) 

PMIN = P(1) —- : 

RMAX = R(1) ; 

RMEIN = R(1) 7 . — 

DO 20 I = 2» 

Me IZ 1. GTe DMIAX DEM £7 (1) 

IF (Z2¢(1)eLTe ZMIN) ZMIN = 2(1) 


I I TT eT ee eee 


oe ee 





MAA 


AANANNDAA 
: i } 
| ie 


686 KOO = KOO + 1 


PAO, 
P(T) 


IF (P(1)eGTe PMAX) PMAX 
IF (P(I)eLTe PMIN) PMIN 
IF (R(I)eGTe RMAX) RMAX = R(T). 
20 IF (R(I)eLTe RMIN) RMIN = R(T) | 
WRITE (69902) ‘ZMAXsZMINsPMAX sPMINsRMAXoRMIN _ 
902 FORMAT(1Xs4HCHE2s6F1004) 
pO 250 I = 1930. 


fou On 


LOCI) = IBLANK i a  —S 
LC(I) = IBLANK ; ede, LN 7 
250 LR(I) = IBLANK 


KHAF =(K+1)/2 
' DO 30 I= 1oK 
JZ =(((Z¢1)-ZMIN)/(ZMAX-ZMIN))*® 290) + le 
JO =(((P(1)-PMIN)/(PMAX-PMIN))* 296) + le 
JR =(CCRCT)-RMIN)/CRMAX-RMIN))® 290) + le —— 
LO (JZ) =ISTAR 
EG (JO) =ISTAR. | 
LR (JR) =ISTAR 
KK = I -— KHAF 
WRITE(65105) KKo (LOCUS) sJJ=1 930) 9 
1KK s(LR( JJ) sJJ=1930) 
105 FORMAT (1XsI4s 1HIs 3OALs1HI s SXsI4selHIs30Als1HI 5X5 
RieslHis30AlsinT )  — i...  ....., 
LO(JZ)= IBLANK 
LC(JO)= IBLANK 
LR(JR)= IBLANK 
30 CONTINUE 


1 
— op ee eo 2 ep 62 om Sp ee oe oe SP ce ce 8 2 2 oe Ss oe ee ew ee 6 6 oe ee ow ow op 6 op 2 Sm eee oh 62> 22S om 6 2D ow SD ow OP SD ow Ces OD ow ew SS ow SD ow ow OR ew | ewe == 


fp ia es es i ee i ee em == om = ee em et ee we ee em ee ee ee ee es —— sa oe ee ee ee ee et ee ee we ee ee ee ce eee 


KK se (LC(JJ)sJJ=1930)5 


ee 


ee me ee ee ee ee ee ee ee ee ee ee ee en ee mm wee me oe es ee at mm mm Ome ome om ome ee me 8 ee et ees 


= Gn ap a =] So SS Se SE SS a Oe as a le ee Se Oe ee ee ee ee ee ee ee ee eg eee eae nee = == 


eigmvO 699 1 = leks 
WRITE(69698) IsZ(I)s P( IT) oRC(T) 
Pee ORMAT(10OXsT593F 2008) oe ee 
699 CONTINUE | ; 
IF (KOO«GTeKOW) GO. TO 400 
GO TO 555 Ba 
400 CALL CLOCK (1) 


> cm OP ce oe ee ce es es i es es es es es ee ee es es es es ee ees ee es ae ee es ee es ees es es ee ee ee ae ee ee ee ee es es ee Se ae ee ee oe ee a SE ee oe eo 


ee ee a meet OOD Ome eam ome ome ome ~ ome om om om ome #- = eee ae ee eh eet ee ne Oe ee ee ee ee es OO Ome ome —- — -- = -- 


_ GROUPING OF C(N) FACTORS (POWER SPECTRUM) FOR FREQUENCY - 
STUDYe POWER SPECTRUM VALUE = A¥®¥A + B*¥B DIVIDED BY —  — 
4.0 AS DEFINED IN THE TEXT 3Y LEEe THIS IS ONLY AND 


ESTIMATE OF THE TRUE SPECTRUM. 





WRITE(69577) | 
577 FORMAT(1H1s37HFREQUENCY DISTRIBUTION OF C(N) VALUES) — 
WRITE (65575) 

575 FORMAT (1H093X91HN 9 6X sGHINT196X s4GHINT2 s6Xs4HINT396Xs 3 =— 
L4SHINT4S 5 6X 9 4HINT5 9 6X 9G4HINT6 9 6X s4HINT 796X sGHINT 89 6X s4HINTS 
26X 9 SHINT10 95X sSHINT1195X sSHINTI12) | 

WRITE (65576) 

576 FORMAT(1HOs9Xs4OHLT 2001 eOO1—7 UL, «OL—s02 M02—=.04 °°  » 
150H-03—-e04 04-205 05-206 006-207 e07T-—-08 e°9 
225HO8--09 009-21 Cicteiinae ———e_o- 





ee ew = ee ee ee ee ee - 


ee ee ee one rd 


eo ee es 


me 213 IF (COEFF. Gis es OB 9950) 


— 


DO 500 N = 
TOT = 0-0 
nO S 
INTI1 
INT? 
INT3 
INT 
INT5 
_ INT6 
INT7 
INT8 
INTS 
_INT10 
INT11 
INT12 0 


1oKKK | 


mt ett | fe me ES SI me 


e 
Oo 


1 | | | | | 


OFM OF OD OD DO 


0 ee a 
0 


th oot 


A oe ooo 


en SR ee AED SF nS eS eS es SD pS GD GE GES weet GD GD GD i GD GD GD De ee OS i ee ee GD Ge ee a ee a ee ee ee 


DO 501 KOO = 15KOW 

GOEFF_= C(KOOsN) = 

TOT = TOT + COEFF 
OTS = TOTSe TE (COERFERCOEFF) 


IF (COEFF.eGTe200099995) GO TO 505 
oy. = INTI + J] 


_ a te te ie ee GED ee S| GNSS RD RE GD SD GD EAD ERS SED cee ees ES ED ND eS eS tes GD GS ED GD eS eS ee ee ee ee ee eee ee ere ee te ewe ee ee ee we ee a 


mm em crm cme mm ee i ee ee et ey ee es ee rm ce ee ec ce ee ee ee es ee ee ee oe 


SseeTO 50] 
BGS IF (COEFF .GTs-00995) GC 
INT2 = 


GO TO 506 
INT2 + 1 
GO TO 501 
506 IF (COEFF .GT++019995) 
_INT3 = INT3 + 1. 
GO TO 501 
507 IF (COEFF +GTe.029995) 


MEGO. TONS07 = 


oe) TO 508 
INTS + Il 
mce- TO. 50) 


608 IF (COEFF eGTe 0039995) GO TO 509 


—a. em ee ee ee em ee ee ee ee ee ee ew ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee aw — a oe 


+ re mn te ee ee — = ee ee is ee et en ee cet eee Se ct ee ce SS ete cee ge cee ce ce ce cee ce ce cee ee ce ce ce re ce ce ce ct te se ee oe as a eee ee 


-— -_— tt ee a ee ee ee ee ee te i ee ee ee ee = ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee et ee ee ee ee ee te ee ee ee ee ee oe 


—-—_ - ee cee eR eh eh pth SS GD Dm re = te te a te a PD ee SD HS mS pS SD SE SD SD SD te i te er ee ee ce cee ee ee ee te ee i ee ee ee ee re ee ee ee ee te ng ee ee ee ee 
' 


we ee eee ee ey ee ee ee — we ee ee ee ee ee te ee ae ee ee ee ee ee ee ee ee ee ee re re ee ie ee ie ee te ee ed a ee ee ee ee es ee ee oe 


—INT6 = INTO. + 1 
mo 70 50) 


(a Rm ert Sues py iene guy i So rms 


510 IF (COEFF eGTee059995) GO TO 511 
INT7 = INTT + 1 
GO TO 501 

_511_IF (COEFF .GT..069995) GO TO 512. 
INT 8 = INT8 + 1 

GO TO.501. 
8512 IF (COEFF eGTe 2079995) 
—INT9 = INT9 +1. | 
GO TO 501 


GO To Siam 


-_— a ee — oe 8 RO ee we 


GO TO 514 


= Oe ee ee te et cr cee me te ee ee me te a ee ee ee ee ee i ee wm es ie ee te ee ee ee ee ee ee ee ee ee ee 


‘ 
oe ee ee oe te — 2 —— io Re ee eet ee eet ee cm tect te cet ete ce ee cei ee es ce cee cee i ee ete i ee ee re ee eye ee ee re ee ee ee ee ee ee OE a ee 


INT10 = INT10 + 1 
feo) 10-50) SS eee 
514 IF (COEFF eGTeel) GO TO 515 

INT] 12]. Niece) een 

GO TO 501 
p25 INT12 = INTI2 +2 
501 CONTINUE : 

AKOW = KOW 

AVE = TOT/AKOW 


-— nee ee ee - te ae ee ee ees cet we | me cent cm cents es cet cee ae te ee ee ee ee ee ee es ee es ee ee ee ote oe oe = 








ee te ee 


TOTS/AKOW 
AVE*AVE 
SORT ( (AVE2- AVE3)) 


AVE2 
AVE3 
SOE V 


LINTOsSINTIOsINTILsINTI2 
502 FORMAT(1HO91T4912110) 
WRITE(69602) _ AVEs SDEV 


602 FORMAT(1HO»s7HMEAN = 9F20¢8s10Xs16HSTANDARD DEVe = 9Ff2008) 


WRITE(69502) NoINT1LoINT2ZsINT39INT4SSINTS SINTOSINT Ts INTB85 


ee ke ~ = sooo ee me creme cr ce crm ee em cm ee en wm ee 0 = ee ee me ee ee et re a ee 


oo 6 oe ee ee ee ee es we es ee ee es es ee es es a es es es oe ec ep Se OD cw oe SE eo a oO 2 a eo 2s GP ap 6 ow SS 2 = oe ow oF a= a= —— 42 ea oe 


= oe oe a a ee OP ee oS 2 Se 2 Se Se ew 22 2 2 eee 2 |e ew ow ce SO ow ce ow ae oe © oe oO 2 2 ee 2 2 2 2 ew 2 ee 2 2 ee Se oO ew ee ew 2 ee ee 2 ew ee 2 OP ow oo 2 ee es ee 


NOM 
m 
=< 
Oo 
oO 
nh 
oO 
A 
Oo 
eG 
Uv 
a. 
G) 
A 
O 
Cc! 
=! 
r=} 
ae 
mm 


500 CONTINUE | 
GALL CLOCK( 2) 


ee ee ee ee ee Fe ee ee ee ee en = = ee ee 


a BE ps Dancearea coneanes sees ceatiroe eens © arin eaten 
CALL EXIT 
END 
SF ae a ares ee oe Sf St CSE ECS 5 SS Sr SS SS SSS ar ISS Se See ce 
SENTRY 
— — ne SS eo eS ee SS a See Sag SS Se Se oe re tn SS ea Sag SS SS SS SS SS a 6S Se Se See ice 
—_ eee ye ye = Se eee Sees 
4 eee ee ery Se eo ro SS Scr oe se eer ee SS oe oF eee 
ee a ee pe SS SSS SoS Sr See Se Sse B® SS hocercore 
= — eee eee aE yy pe 
é 
— Se a eae >. ae a ee 
ee an ete i . e Sse SSS Sp SS ara SS 5 SS SS SS SSeS Sr SSeS oe ee ee 
= — = ae ee a tea cone ra ee ——————————EE a = ees ae oe 
= — = = S| = Ss See Se Se a ate ote ee eee oe a om te ee ST ST QS SY Te ER ee te SE CD GS GD ep cm gee en 6 ee ee cm we em SS 
= Se ne ewe or a pe ee — = see oo SSS ES eS a eS SS SS SS SS St So SSS Sr SS Sa SSS SS eS Be Se Sree 
4 
— Es crsennennsges oan a a a ee 
Sos See = re Se Se See se SS ccm cen Re. erst SS or Se ear Sea Sea SS Sa tt SS SS SS 


— a ens oem em ete ame 
— _—-+ Se Ra a eee ee - me ee pe es cm ce a ee ee ee cn ee me ee ee ee ee a se Se -- 
é 

= -_— _ —— ’ ep me ec em te ec ce ee ee es + — ee Sk ce A mS ms ee gm ms ms ee es cm es ct mms em ee See cee es ee gee ee ce ee ee ce ee ee ee ee ee ee ee re re 
ee See a eee ST 

— et a 21° ee ome ee ee ee ee ee ee ee ee em ewe + os moe = 

7 
— = ~_ — oe oe _ — ae SS ee oe wee ee OS ee — 
——: es senreren erte = me he en em ee eee eee 








i * 
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SINGLE FOURIER 3ERIE3 PROGRAM 


Le 


| READ: TITLE | 


READt K,EL 


KOW = No. of 3Jamples 





| 


KKK = No. of Harmonice 


DELX = PIE/(K-1)/2 


Z3UM = O- 






555 Frou 
End of 
Programe 


SUMOS = 0. 


z(1) = 0. 





X(III) = X(III-1)+DELX 
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(4) 


( READI AZ¢ | 


, 





= se 
ui 2(dd)= 1. | 





i ’ 


ge Fe es ell 


} 
| 
| 
E 






| SUME =(2(1) + 2(K))* 4 
| AZBRO = (2./ANK) * (730M + 3UME) | 












SUM3IN = 0. 






SUMOO3 = 0. 






SOME = 0. 






FN = N 


3UM =< 0. 






SUMOO3 = 3UMOOS + 2(I) * 003 
((FN * PIB * X(1))/8L) 


QUMSIN = SUMSIN + 2(1) * SIN 
((PN * PIE * X(1))/6L) 
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— —— ae 


| 30ME = ((2(4) = 2(K))/2 * UOS(FN °¢ srs) | 
nego ar 





| 
At 


A(N) = one ° ( sumoos + 30NE) ; 
H(N) = 2./AlK * (SUNSIN) 
[0 oc = aN e aK) + a(n) ° “any | 7 


1o(Koo, N) => 00/4. 


ao eee wn 2 ane | 


WRITE: O(KOO,N) 











NMAX = K-1/2 
’ (1 Fi) 1 
| 
: SUMC3 = A’BRO/2. { 
| 
(a > ae 
{ 
a 
| 
FN = N | : 


SUMO3S = 3UMO3 + A(N) ® 
COS((FN * PIE © x(1))/ 
EL) + B(N) * SIM (FN * 
PIB * X(1))/8L) 


Ones 
| 


P(I) = 3UMOS 


x A(I) = P(1) = 2(1) 


33 


— 


PLOTTING aout] 





CALL exIt | 


—— a 
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——— an 


SINGLE FOURIER SERIES rLOTTING ROUTINE 







From Main Program 


a ee 


ISTAR = 13702990896 







IEYE = 27661634608 
IBLANK = -17997958192 
IDASH = -818089008 





) tosHcx) = IDASH 


i | 
\ on 


peaites IDASH LINES 


bn 


an, 


~-2MAX = 2(2) 
ZMIN = 2(1) 
PMAX = P(1) 
PMIN = P(1) 
RMAX = R(1) 
RMIN = R(1) ; 


ie TEL 


<> 2 fos 
Se 


> 


i a = P(I 
) 
’ 









AMIN = R(T) 


WRITE: 

ZAAR, ZMIN, 
PMAX, PMIN, 
RMAX, RMIN 


LO(I) = IBLANK 
LO(I) = IBLANK 


LR(I) = IBLANK 


| 
' 
{ 
ee cae er na onic) 


] 
| 
L 


KHAF = K +1/2 


Capa 


JZ @ (((2(1)-241N)/ 
(ZMAX-2MIN)) * 29.) +1. 








JO = (((P(1)-PMIN)/ 
(PMAX-PHIN)} * 29.)41. 











JR = (((R(1)-RMIN)/ 
RNAX-RMIN}) * 29.) + l. 









LO(JZ) = ISTAR 
LO(JO) » ISTAR 


LR(JR) = ISTAR 
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+ 











KK 


I=-KHAP | 


YRITE: KK, 
LO(JJ), KK, 
LO(JJ), KK 
LAC UJ 


















ITtLANK 





LO(JO) IBLANK 


InLANK 


—_  — 


RETURN TO HAIN OUTINE 


-----@ 


— = ee 
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COMPUTER PROGRAM TO DETERMINE COEFFICIENTS 
FOR DOUBLE FOURIER SERIES (TWO DIMENSION) 
APPROXIMATION OF A SURFACE AREA 


Purpose 


The purpose of this program is to approximate a surface 
area, f(x,y), using a double Fourier series, and to present 
the Fourier coefficients which are a characterization of the 


surface. 


Language 
Fortran IV (IBM 7040 Computer) 


SM OCi tem mre U1 Ona iy 


Variable 5/A* iy Ores Description 
X A I Independent variable data 
point along x-axis. 
nd A | I Independent variable data 
point along y-axis. 
Z A I Dependent variable data point. 
EL S il One-half of fundamental per- 


iod along the x-axis. 


mC S lise One-half of fundamental per- 
iod along the y-axis. 


KK S I Number of data points in 
x-direction. : 


* S--Single Variable; A--Array of Variables 


** T--Input Variable; O--Output Variable 
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Variable  S/A iO Description 
LL S Ne i Number of data points in 
| y-direction. 
MM AX Maximum number of harmonics 
and S I to which Fouricr series sismcom 
NMAX be expanded. Should equal 
| 7 one another. 
M 5 Oa Harmonic index. 
N Ss O Harmonic index. 
Le S If Mathematical expression input 


as being equal to 3.14159265. 
Fourier coefficient. 
Fourier coefficient. 


Fourier coefficient. 


> Fe FP PLP 
So oes © 


Fourier coefficient. 


Se ee, © San 


Program Routine 
This program utilizes data points equally spaced upon 
a grid, placed upon the surface to be analyzed. An-approxi-. 
mation of the surface is created by expanding a double 
Fourier series of the type discussed in Appendix B. The 
Fourier coefficients are presented as characterizations of 


the surface. 


Input Generator 
Immediately following the listing of the double Fourier 
series program, is a listing of a simple data generator for 
the program. This program uses a random set of coefficients, 


and Equation (2) of Appendix B, to create, in the form of 
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punched output, data points simulating a surface area, which 


are the appropriate input for the double Fourier program. 





——— ae a et at an 
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PROGRAM 2D FOURI _ AUGUST 1965 
DOUBLE FOURIER SERIES (TWO DIMENSION) PROGRAM 
PROGRAMMED BY WD ALDENDERFER _ ; cm 
DIMENSION X(25)9Y(31)92(25931)9A(595) © <i 
DIMENSION EN AS) | UGA BER I 2a 0 | en lei - 
CALL CLOCK(0) 
WRITE(69600) _ : a ol” 
600 FORMAT(1H1»19H2D FOURIER ANALYSIS) 
WRITE (65623) | | - le ee «oe 
623 FORMAT(1HO0920Xs26HDISPLAY OF CALCULATED 2-D +4 
LS4OHCOEFFICIENTS WHICH CHARACTERIZE FUNCTION) — ie 
WRITE (65654) 
654 FORMAT (1HO94Xs1HMs4X91HNs 7X s6HA(M9N) 9 14X 96HB (MON) 9 7 
114X s6HC(M9N) 914X96HD( MN) ) _——— tt 
READ(5592) EL sHsKK os LL 
2 FORMAT(2F2008s 215) ia a ee eee ~ ak oe 
g Et = FUNDAMENTAL LENGTH IN THE X DIRECTION 
C H = FUNDAMENTAL \LCENGTH IN THE Y DIRECTION © - es 
‘ay KK = NOs OF DATA POINTS IN X DIRECTIONS 
C- LL = NOe OF DATA POINTS IN Y DIRECTION’ a 
-READ(551) ((X(TI9Y(S)9Z( 195) 9 J=LoLL) 9131 sKK) 
RIEORMATI3F20¢8) em eee a 
C X = INDEPENDENT VARTARLE. 
Memes ¥Y = INDEPENDENT VARIABLEs = 0 
C Z = DEPENDENT VARIABLE. 
| MMAX = 5 = i S - 
NMAX = 5 
C  NMAX AND MMAX ARE THE NUMBER Soe ARNON CS SO. a me i 
SERIES IS TOsBEwEXPANDE Dj. semen cern ee 
BIE = 3.147592765 Ge lee 
AKK = KK-1 es ee eee ee 
ALL = Lt-1 F 
MKK = AKK Eee oe es i ae 
eis = AL eee 00 ae 
AMK = AKK#ALL 
DO 22 M = l»eMMAX  —. oO age — a ia nna. 
FM = M—-} —_ | 
DO 22 N = leNMAX i a rr a 
FN = N-1 Se mo ae 
ABLE = Oe 
ABAKE = Oc 2 ee a a 
ACHAR = 0.0 
AD EU lueten oO _ a 7 axe en = 
RCHAR = 0209 
BDFLT = 0.20 7 | ee 
CRAKE = 020 — eS 
CDELT = 0.0 
Sorel = Oso ; _ 
PIM = PIE #* FM 
PIN = PIE#FN as “a 











— — ee oe ee ee ee 
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C ee. : : 
rc CALCULATE A(090) 
C 
| ABLE 2° (0 CC(Z01 Sm 2°01) ee ee, CO 1) 20K 
me ~—Ss1/2.0)} * COS(PIM) * COS(PIN}) a= _— 
DO 20 I = 29MKK 
’ DO 20 J = 2oeMLL _* «+ See a 
ARGX = PIM * (X(I)/EL) 
ARGY = PIN * (Y(J)/H) nc Deen 7 
CX = COS(ARGX) 
CY = COS(ARGY) ——s ne ae — 
SX = SIN(ARGX) 
SY = SIN(ARGY) | oe eee ; 
IF([1eGTe2) GO TO 400 
C 
C CALCULATE COEFFICIENTS BY PARTSe 
C 
ABAKE = ABAKE 4#+((Z (194) + Z(KK9J))/200) ¥CX ¥CY 
s GBAKEe= CBAK Em ilileZ t) oJ eatnZ (Keanna). iG Xl SNe ; 
BPOOUIF(S.-GTe2) GO"TO 401 
ACHAR =: ACHAR+ ((Z(Io1) + ZC TsLL))/200) *CX¥CY sss 
BCHAR = BCHAR + ((Z(19o1) + Z¢(TsLL))/260)¥*SX¥CY. 
401 ADELT = ADELT + Z(IsJ) *CX*CY | 7 Ls 
BDELT = BDELT + Z( I oJ) *SX#¥CY 
POEL T = )GDELT + Zito te OK SY ee 
ODELT = DODELT + Z(IeJ) * SX * SY 
20 CONTINUE. i __._. Se 
A(MeN) = (ABLE+ ABAKE+ ACHAR + ADELT) * (4e0/(AKK*ALL ) ) 
— B(MsN) = (BCHAR + BDELT) *(4¢0/AMK) - : 
C(MsN) = (CBAKE + CDELT) * (46¢0/AMK) 
mOeMiMeN) = (ODELT Lm Se07ZAMe = UCC (‘(‘(‘#’(UUCOCOOUUUU.U Uo 
| 22 CONTINUE 
CC eee eae: eee. 
C WRITE VALUES OF COEFFICIENTS.« 
C 
DO 25 M = 1»+MMAX 
— MM = M=-] en eee . 1s 
DO 25 N = J »NMAX 
IN NS a i no 
WRITE(6927) MMs NNoA(Ms N) o+B(MoN) oC(MsN) sD(MoN) 
27 FORMAT (1H0O921594F 2004) mi. ; : ; 
25 CONTINUE | ; 
| meat  CLOCK()) ee MS. | ae ee, ne 
CALL EXIT 
4 (2) ee ee ae eee. ee... oe 
SENTRY 


ee ee ee ee ee ee ee ee ee ee 





SM) OS SS ie i a ie 
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PROGRAM 2D GENERATOR 
PROGRAM TO GENERATE DEPENDENT VARIABLE Z AS A FUNCTION 
OF TWO INDEPENDENT VARIABLES X AND Ye PROGRAM 
USES DOUBLE FOURIER EQUATION AND GIVEN COEFFICIENTSe 
COEFFICIENTS USED IN THIS PARTICULAR PROGRAM ARE 
TAKEN FROM PRESTON-HARBAUGH PAPERe 
PROGRAMMED SEPTEMBER 1965. 


DIMENSION A(595)9B(595) 9C(595) 9D(595) 97(25931) 
CALL CLOCK iM eee _ 
DO 100 N=195 
DO 100 M2195 
FN=N 
FM=M 
A(MsN) 
B(MsN)- 
C(MoeN) 
D(MsN) 
100 CONTINU 
DO 101 
C(Mosl) 
B( 15M) 
D( 1M) 
101 D(Ms1) 
m= 3414 
SUMZ = 0.0 
AH = ia 
AL = l2e 
WRITE(659984) 
Bee FORMAT(1H1s29HTNPUT VALUES OF COEFRICIENTS*) === = 
WRITE (69982) | 
982 FORMAT(1H095X981HM §N ca LM oN) B(M9N) 
l C(MoN) “D(MoN)) 
DO 300 N = 195. 
DO 300 M = 1595 
WRITE (659981)M9NsA(MoN) sB(M9N) 9C(MSN) 9D (MN) 
PS lMPORMAT (2Xe2159GFcUeC) we We 
300 CONTINUE 
PIL=PI/AL 
PITH=PI/AH 
XC = -l2e 
YC = -156 
INC = Ted 
WANG = 1.0 
WRITE (69983) 
983 FORMAT (1H1 911X924HX(1) le) Jap 


ee i a a ee ea on — 


1?2H ZC Yoo) ie J) 


, 0 eee comm cms me mee 2 ee es cm ee ee ee es ee es . ee ee ee = oe 


ee ce meee ee ey ce ce ee cre meg ee eg ee ee ee — - os — ee oem oe -_—_— = 


“Wed + eO2kEN + « 1] L*EM) E00. | 
(o2 + eO4*FN + 208*FM)*100— ea caiaiicinss ee 


(ol + eO6*#FN + «O5*FM)¥1000 
(208 + 0 O5S*FN + 002*FM) #100. 


we ae ee om oes et — ome o 8 ane oe : oe ome ome eg ee oe 


nod nf on ST th th 


—— os = =e 








a — oo a a a 


ene a eS OO EE © ee om ae: a a SE a RT a a TR: ah 


a oe ee ets ss ee = oe es oe © + eee a ma ee ee os meme cm cm} arm are any eee gm ey my rm ee eg eg ae ey weg es weg og lh wee cee cm oe oy oy oe ee ow | ne cme ree cee es ee ee ree ee ee es ee a we es es en oe es 


oe ee te oe ee ee eee 


ee ee ee ~~ _—_ oc > ame ams 


PINH = PIH*FNT 





Ee os oi 


DO 200' M=195 
FMT=M—1 
PIML = PIL*FMT 
IF (N=1)14091203140 
mee 1F(M—-1)15099305150 aaa 
130 AMDA=e25 ; 
GO TO 170 ae ree a I. ON ck tee ae. 
140 IF (M-1)16091505160°. 
150 AMDA=.50 
GO TO 170 
160 AMDA=1-00 .- 
170 ARGX = PIML*X 
_ARGY = PINH#Y 
CX = COS(ARGX) 
= COS(ARGY) — . 
SX = SIN(ARGX) 
= SIN(ARGY) —_. 
A(MsN) ¥CX¥CY 
BI MsN)*#SX*¥CY 
C(MoN)*¥CX*SY 
D (Ms: NHS XOES Ne ee 
TS = AMDA*(ACC + BSC + CCS + DSS) - 
q Pe = ZC + TS 
200° CONT INUE 
_ WRITE(6s980)° Xs¥9ZGoil oJ m0) Pee eens ee 
980 FORMAT (1X%93F20¢83215) . 
a WRITE(79950) XsYoZC ———s—“—ssssSsai‘(‘( ai ——————— ts 
950 FORMAT (3F 208) \ 
a LE lh MO ee 
Oe 
YC + YINC 


et cy te tee SA ey Se ey ie cee GE cee cere em 6 me cm mm ce cere cme cee cee ce cee cee em cee ee ce ee ce ee ee ce cee ee ee ee ee ee et ei re ee ee i ee ee ee ee ee ee ee we ee ei ee ee i = 


ame e-< soe ae — 2 - i a a Sk ee Ge GS GE GS AE ie ee ee ee ew ee ee —_ —— ae ee eee 


eae ieinn _—— wae — ey me ey me ce ee ey ey me ee ee 2 ee ee ee ee we ee we a ee re ey ee tw ee a ee 


Se) 
Ca) 
CY 
| | 





5 
© 
wm 
< 
CY 
“th oth ou 


Bewcmxc = XC + XINC =~... h...ULULULULUDLULULULULLULULUL oe 
T(lol) =(TE191) + T1931) + T2551). + T(25531))/4e0 
Mee (1531) = 020° ~— = 
T(2581) = 040 
T(25931) = 000 | 
DO 100s¥ = 93 
_ Mle S)= (Talis 
T(259J) = Oc 0 
700 CONTINUE | 
Do 7OY | = De 
meee C19 1)-8 (TCI 91) + Tlls31) 1/2000 US 
T(I1o31) = 040 
__ 10] CONTINUE 
DO 702 I = 1925 
0O 702 J = 1931 
SUMZ ="SUMZ + Tlels J) | 
T02 CONTINUE 
AO0O = (460*SUMZ)/7206 
J WRITE(699999) SUMZs AOO | 
9999 FORMAT(1H1s18HSUM OF Z VALUES 
CALL CLOCK (1) 
490 CALL EXIT 
: END. 
.. BENTRY 


_— — -_—— = —_ 


oe ete ee cee et ee ee ee ee -e _ ee rw me me eg ey ee ee ee en ce ee es ce ee wr crn ee we ss mee me sw Ss 





ee ee ee eee ee ee es ee 


“9F200895Xs9HA(090) = 9F 2008) 


es 6 ee 





ee — 


SO ee ee ee es es es oe = oo ee i ee ee ee ee ee ee ee te ee ee ee ey ey ee ee we ee ee ee + oe a=: — _ _— 





DOUKLE FU RISER SEATHS PROIRAM 





[ vanes PI fe 
\ and HEADINGS 


~ 


’ aN j 


( R5AD8 BL, H, AK, tL | 


ReaD: x(I), ¥(J), Z¢ 1,9) | 






N4AX = NO. HAAMOD ICS 


NO. HARMONIO3 





| PIE = 5.14159265 





ADELT = 0. 
| ABAKE = 0. 
~—ACHAR = 0. 

— -BOHAR = 0. 
| BDELT < 0. | 
; OBAKE = 0. | 
| QDELT = 0. 

| 


' DDELT = 0. 


oe) 













PIM = PIES * FH 


PIN = PIS * Fh 


| 


(1,1) + 2(KK,1))/ 
1,LL)+ 2(KK,LL))/ 
¢ 903 (PIM) * 003 
























OX = CO3( ARIX) 
CY = 003( ARGY) 
3X = SINC AROX) 


39¥ = 3IN( ARGY) 


ABAKB AFAKBE ((2(1,0) 


OBAKS OBAKE ((2(I,J) 
Z(KK ,J))/2.0) * Ox * 8 


ACHAR AQHAR ((2(I,1) 


BCHAR BOHAR ((2(I,1) 


ARIX = PIM * X(1I)/BL 
ARGY = PIN * X(I)/H 


Z(KK,J3))/2.0) * Ox * OY 


2(I,LL))/2.0) * Ox * oY 


2(1,LL))/2.0) * 8x ® OY 


2.0) 
2.0))/ 
(PIR) 









Y 


96 








ADSLT 


it 


“OELT 


i} 


ONELT = 


tt 


DOSLT 


ADELT +701 ean eee 
HDELT +t ZC 1,5) * 3X * OF 
ODELT + 7.(1 J nS oF 


DOSLT+2Z(1.0) * 4% * SY 





4(M,N) = 









(AHLE + ABAKBE + AGHAR + ANSLT) 
© (4./aKK * ALL) 


B(M,N) = (ROHAR + HNELT) * (4./ AUK) 


— 2 —<——=- <b embans 


WRITS: M4, NN, A(M,N), 
B(M,N), O(M,N), D(M.N)], 


ZG 
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COMPUTER PROGRAM TO DETERMINE COVARIANCE 
‘ FUNCTION AND POWER SPECTRUM 


Purpose 
The purpose of this program is to determine the 
covariance function and the Sone EE ECan of a one- 
dimensional function, f(x). The secondary purpose is to 
plot the power spectrum and the confidence bands of the 


Pmectrum, CO facilitate analysis of the tfunetion. 


Language 
Fortran IV (IBM 7040 Computer) 


SyVMloOLLemune FlOonary 


Variable S/A* 1/0** Description 
a A I Value of the function, at each 


data point. Read in as fixed 
CON ENB WINS 6 


| ~ 

if 5 -- series index. 

K S -- Lag Index. 

Go A O Covariance estimate. 

L A Osa Raw estimate of the power 
Spectrum. 

U A O smoothed estimate of the 


power spectrum. 


* S--Single Variable; A--Array of Variables 


** T--Input Variable: O--Output Variables 





29 
Variable o/A I/O Description 

N S ik Number of data points. 

M S if Maximum lag. (Less than N/3) 

EDF A it Equivalent degrees of freedom. 

TL95 A I Table of confidence limits, 
lower 95 per cent. @ 

TLIO A | I Table of confidence limits, 
lower 9O per cent. 

TU95 A it Table of confidence limits, 

| upper 95 per cent. 

TU9O A If Table of confidence limits, 
upper 90 per cent. 

KOW = I Number of samples. 

KOUNT S -- sample index. 

LA A I Storage for data identifi- 
cation. 

LAG S O Degree of offset by which 
the function is correlated 
with itself, to determine 
covariance estimate. 

UMAX 5 -- Maximum value of smoothed 
power spectrum in each fre- 
quency band. 

UMIN i) -- Minimum value of power 

. Spectrum. 

SPANA S -- 90 per cent confidence band 
of smoothed power spectrum. 

SPANB 5 -- 95 per cent confidence band. 


Program Routine 
Thi's program utilizes data points which represent the 
value of the function to be analyzed at equally spaced 


* 


intervals along the function. Using a numerical integration 





100 


technique, the trapezoidal rule, the values of the 

covariance estimate are determined according to Equation (3.25) 
of this paper. A rough estimate of the power spectrum is 
computed using Equation(3.26) am smoothed according to 

Equation (3.27). The table of confidence limits presented by 
Granger’ is read into the program and the ninety and 
ninety-five per cent confidence bands are determined. 

Ms smoothed power spectrum with these confidence bands is 


displayed graphically. 





| 


I 


' 


AN NNNANDNNANNANNNANNINNANIANNANININANNANINNANAINONAAN 


DIMENSION X(3200) 9C(301) ok (302) 9U( 301) sEDF(13)sTL9O5(13)5 | 
1TU90(13) 9TU95(13) sLA(13) »JOW1 (53) »JOW2 (53) 9TL90(13 ) »MX (3200) 
“REAL L 


-_—— a aw ee awe OP a= ap oe ce cee om oe coe SF oe ee ce oe Se oe ce ees es es ee es ee ee ee es ee ee ee ee ee ee ee ee a eee eee ee ee eee Ss ee ee a ee ee a ee oe 


PROGRAM COVAR FOR COMPUTING COVARIANCE FUNCTION 
AND POWER SPECTRUMe PROGRAMMED SEPT. 1965 BY 
FW PRESTON AND WD ALDENDERFERe _ 





NOMENCLATURE 

= LAG INNFX =. tne le... ee 
I = SERIES INDEX 
J = WEIGHTING FACTOR INDEX 2 __. A 
C(K) = COVARIANCE ESTIMATE 
L(J) = RAW ESTIMATE OF POWER SPECTRUM ers ee 
U(J) = SMOOTHED ESTIMATE OF POWER SPECTRUM 
M = MAXIMUM LAG (LESS THAN N/3) ss ——i‘(‘(‘( 
N = NUMBER OF DATA POINTS 
LA(I) = STORAGE FOR DATA IDENTIFICATION ————C Ns 
EDF(I) = EQUIVALENT DEGREES OF FREEDOM 
TLO5(1)= TABLE OF CONFIDENC LIMITS - LOWER= 95PERCENT _ - 
TLOO(T)= — LOWER- 9OPERCENT 
TUSO(T)= em UPPER 90 PERCENT) 
TUSS(T)= - UPPER 95 PERCENT 


KOW IS THE NUMBER OF SAMPLES _ 


ee a a ey ey ec ee es ee ee ee ee ee ee ee ee ee ee ee ee ee —— a i a a  — a 


ee ee em ee ee ee ee ee ee ee oe ee ee om oe oe eT TT a Ne ee ee ee 


ern CLOCK....f.0) 
KOW = 2 a 
KOUNT = 1 : 
READ (59110) (EDF(T) sTLO5(1) s TLIOO(I) s TUIO(T) » TU9S5(T) 9 T=1913) © 


ee ee ee ee ee ey ee ey oe I I eT me ey ees ee ee tes ee ey tes tes tes tes tes es ey es gies tee its ets gies ee ee ee eee ee tee po ee Oe Oe ey ey ey ey ey ey ey ee ee ee ee ee oe eee oe ee ee ee ee 


110 FORMAT (F501 94F5e2) 


ee ce ce ce ee 0 ee ee ee ey ee tte ee ey ey ey ei ey my ees etn me ee em tems ies ete tee ee re ee ee ee ies tee tee te oe ee ee ee ee 2: i ee i ee ee ee ee ee es ee eg ee ee ee oe oe oe 


400 READ (5599) (LA(I)s1T=1913) 
99 FORMAT (13A6) 


eco ee 


me mm ct mm mm me we we rs th | a ee ee ee we i es ee os cc i we es we em se ee ee = ee 


READ (53100) NoM 


100 FORMAT (215) 


700 CONTINUE 


et 


101 FORMAT (8011) © 


—_ — —— mmm i I RR a ae NR 


READ (59101) (MX(T)»9T=19N) 
BO 70081 I= -IsN  ~ et ae I ee 


X( 11) MX( IT) 


at ee cy cme em ee ee ee me ee ee re ey a ee eee ee ee ee ee ee ee ee ee Oe ey et ee ey ee ee ee ey ee ee ee ee em 8 om om am ss ay oe es ot 


Ome) On 1 = LsN 
MeO TY = O% Jllslls10 9% 


mA 1) = -1.0 - 


10 CONTINUE 


FM =M 


Se ce es DS ees es es es ee ee ee ee 9 De Oe ee ee 4 85> et ye i em es es ee ee ee ee ee ee oe 


COMPUTING COVARIANCE ESTIMATE C(K) 


en mm rem tt mt ee ee Ee I Re om ee ER SRE 


MP] = M+] 

FN =N 

J = 0 ~— 
AK = 2e0*FN/FM 

DO FLO Mise Masel? Caled, 

J = J+l 
IF(AKeLTeEDF(T)) GO TO 320 


—— — my Ge ey ey ey ay ey _ — yy oe ee ee ee 


910 CONTINUE 





aa =) JS =" 
CHG = (AK=EDFdg)eeAGEDF Gos =EOr i): } 





“T9OL = TLOO(J) + VTEIO( Oe Ce CO : 
TO5L = TLOS(J) + (TLOS(J4+1)—-TL9O5(J) ¥CHG 
T95U = TU95(J) + (TU95(J+1)=TU95 (J) )*CHG ae : : 
T9OU = TUSO(IS) + (TUID(J+1)-TU9IN (J) ) *CHG 
DO 210 KKezeles MON eee eee 
K = KK] 
Sl= 0.0 Be ee bs we eee 
NMK =N=K 
FNMK = NMK _. 
KP] = K+l 
DO 205 T=leNMK - re os ; —— 
IPK =I+K | . | 
FO5=S1 = S). + X(P)EXCTPK) eee 
S2 = 0-0 é 
DO 206 I=KPI1.6N ee UL eee J. 2 
206 $2 = S2 + X(T) 
_ 53 = 0.0 ; = es. A __. a 
DO 207 T=1 o»NMK 
ZOg4S3 = $3 + XI) sj _ . Ses >... 22°50 A fee 
210 C(KK) = (S1 — (S2*S3)/FNMK) /FNMK 


“COMPUTATION OF RAW ESTIMATE OF, POWER SPECTRUM L(J) OO 


re 





MN1= M=1 
BO230e3) ="15NPIe Fl t—(“i‘“‘“<“‘<‘i‘ia 
J= JJ-1 
FosJ 
$1 =0.0 
DO 220 K=1l»MN1_ 
FK=K 

220°S1 = S1 + C(K#1)*COS(Z3e1G415926*FK*FI/FM) 


230 L(JIt1) =(260*S1 + C(1) + C(M41)*COS(3614159268FJ) ) 4602831952 


L(¢21)2L 03) 
L(M+3) = L(M4+1) 


mm a ne cee cre te ee oe ee oe re ne ee i a ee ie a i 





eee 


ac R ee 
C COMPUTATION OF SMOOTHED POWER SPECTRUM Sik glee 
a Se ee een 
MP2 = M+ 2 
_ e240 ee ZwR2 gee as. 
eGo Ut) = .25*L(J-1) + .50*L( J) + 625*t(J+1)} 
oo a ae ee a o 
g WRITE HEADINGS ; 
PG - lt. i GOI a 2 


— WRITE(69103) 
137X%s24HFOR DISCRETE TIME SERIFS/31X%s—t—CSsstSttS 
237HUSING TUKEY HANNING WEIGHTING FACTORS///) 


WRITE (6910G) (LACT) 9121913) | ie eee 
104 FORMAT(1X913A6///) ; ee le ee 
WRITE (69106) 
= hy 7 = . Pn ee Bes 
C WRITE POWER SPECTRUM AND TOLERANCE ESTIMATES 
C 


RA EET, RS | A LT TT A A nt mn te re on Ey NN = —_ me 





CV OVO’ 


250 
FOS 


T64 


207 


106 


NADA 


DO 280 T=1953 


280 


108 


290 


300 


— BLOHCOVAR TANCE 95X 9 3HRAWs 9X » SHLOWER 9 19X 9 SHUPPER » 7X9 


LO3 


DO 2 > Ure =2 sme 


LAG = I= 2 

U9O0L = T9YOL*#U(T) E 
U95L = TO5SL#UCT) 

U90U = TINUFU(T) ; 
U95U = TI5U*#U(T) 


WRITE(65105) LAGsC( 11) sL (1) sU9OL UIT) sU9SUSUISL sU(T) 9U9OU 
FORMAT (1X91398(2X9F1004) ) 





PLOTTING ROUTINE 


ISTAR =-13702990896 ~ 

IMINUS = 29809118256 en 

IPLUS = 17997958192 7 
IBLANK = -17997958192 

IT = 27661634608 ; _ : 


WRITE (69103) 

WRITE (69104) (LA( I) sT=1913) 

WRITE (69107) 

FORMAT(19Xs23HSMOOTHED POWER SPECTRUM»34X59 | 
123HSMOOTHED POWER SPECTRUM /15X>¢ 

230HWITH 90 PCTe CONFIDENCE LIMITS s27X9 © 
330HWITH 95 PCTe CONFIDENCE LIMITS/6X%-o 
451HO0O l 2 3 4. is, 6 vo) OF 2 oO 1 Os 
56X »51HO 1 2 6 4 5 6 i 8 9 0) 
FORMAT(LHO939X5 
LS5SOHSMOOTHED SPECTRUM WITH INDICATED CONFIDENCE LIMITS/ 
234X97THIO PCTesl7X97THIO PCT ¢s5Xs7TH9O5 PCT esl 7X9 7HO5S PCT e/8Xs 


oe et ee ee a ee 





cee ee ae ee 


56X sSHLIMIT 6X9 7HAVERAGE 9 6X sSHLIMI Ts 7X9 5SHLIMIT 96X99 7HAVERAGE 9 
66X »SHLIMIT) a << nme aay ee 


PLOTTING OF SPECTRUM 


2s = - — ae < -- o—. * See Te Te de 


JOW1(1)=IMINUS 

JOW2(I1})=IMINUS  — a - 

WRITE(65108) (JOW1(1)591=1953) 9(JOW2(1) 9121953) 

FORMAT (5X953A195Xe53AL) (1°. 00 nner a 
DO 290 1=2552 
JOW1(1)=IBLANK 
JOW2 (1) =IBLANK 
JOW1(1)=T71 | 
JOW1(53).21 
JOw2(1) =! 
JOW2(53)=1 
“UMAX = U(2 
UMIN = U(2 
DO 300 I = 3sMP2 

IF (UCT) eGTeUMAX) UMAX=U(T) 
IF (UCT) eLTeUMIN) UMIN=U(T) 


ee ee ee ee ee ee me Ore Ore Ore > <> ne me ee ee oe ee OE GD Bam ome re ee gee ee Oem Ces ee Or foe OD eee eum 2 orm cee Od foe ot etm eum ee ee em a ee 


——s - eR ee RT I Se hte 


CONT INUE 
SPANA = ALOGIO(UMAX*T90U) ~- ALOG1O(UMIN#*T9OL ) 
SPANB = ALOG1O(UMAX*T95U) = ALOG1O(UMIN*T95L) 


ee — 


Sere a<metr a 


> ay a ee ee Sree wee ee ee = 


See ee a ee = 9 age Sabo oft cit cin aA cic ahh nth GOD ea tee ate! atch nish eum anh =b- ah Ghar Chm as Gun que Gun qu am qus aun ame umm sb Gam ams ats oan an am ok am ae amam os as < 


— i > <> oa eee ee 


mm mee ee em ee 


~~ 62 





ZA = 
ZB = 


DO 3 


IA1=((ALOG1IO(U(T)*T9OL) - ZA)/SPANA)#50. + 200. 














ALOGLO(CUMIN*T9OL ) 
ALOG1LO(UMIN#T9S5L ) 
10 I =2 ¢MeZ 


py etn RR NR pg 


TA2=((ALOG1O(U(1))-ZA)/SPANA)I*®50¢ + 20 


IA3=((ALOG1O(U(I)*T9IOU)— ZA)/SPANA) *506 + 260 
— JOWL( TAL) = IMINUS © 
JOW1(1A2) =ISTAR 


ee ee ee ee ee eee ee ee 


ee es ee ee 2 EE SS ep 


JOW1(1IA3) =IPLUS 


IB1=((ALOGIO(U(I)*T95L) - ZR)/SPANB)*500e+ 240 


IB2 


IB 3=( (ALOG1O(U(1I)*T95U)-Z2B)/SPANB)I*50¢ + 200 
JOw2(1B1) = IMINUS 
JOW2(1B2) = ISTAR 

JOW2(1B3) = 


LAG 
WRITE (69109) LAGs (JOW] (J) sJ=1953) sLAGs (JOW2 (J) »J=1 953) 


109 FORMAT (1X91391Xs53A191X91391X953A1) 
JOWL(TAL) 
JOW1(1A2) 
JOW1(1A3) 
JOW2 (181) 
JOW2 (182) 
JOw2 (1B3) 

310 CONTINUE 


NAO 








= ((ALOG1O(U(T))-ZB)/SPANA)*50¢ + 240 


—— ee ee cy ey ee ee ee ee eee ew oe 


_—_ aa a ee ee em me ee _ 2 eS — ~ ee a Ne 


IPLUS 











—— « -_—-. eh ey ny et eh — — = om ~— — m= « — - cme cime cteme setie ae 





=el 2 


— oe — - ~~ ro Pa me A ey ey ee 


-IBLANK 
IRLANK 
IBLANK 
IRLANK | 
IBLANK 
IBLANK 


| | | 


ee ek ety ee met ae ten mete em enh em my etm ey em ee tts ten eth th cme mm emmy ty im es ete ee ee ie i ce cee GS ate GS ete ie cee we ee ee ED Cy SD Sey Cy ED GS OD ee ey ee ee ee ee ee ee ee ee ee ee ee 


167 IF (KOUNTeEQeKOW) GO TO 765 


ee 


= ae ce tect cee te et ee ct cee ee ee ee ee ee ee oe ee ee i et ee ee ee ee ee ee ee et es en ee et ee ee ep ep ee ep et ety tp te ee ee ee ee ee ee ee ee ee ee ee = 


CALL CLOCK (2) 


KOUNT = KOUNT + 1 


fe eee ee 


~ GO TO 400 


mroswecAtl EXIT 
END 


- 


l2e 


206 
30. 
50 « 

— 1006. 


— _—— a 


<eCOW 


027 


0 34 


039 


0 43 
2 48 
054 


062 


«69 


— re 


me me tm ee ee ee me 8 ee ete ete em ee ce ee ety ee ee ee ee ee ee ee ee ee cee ee ee i ee me =o 


ws 2ieBiOr2 «99 
720 2608 2860 
612.6 nein 26 Sat 
32 1685 2621 
1 AN a ee ees 9) 
044 12668 1494 
249 1660 16683 5) “SU R00 00t0t0 00000 

ace WANG) 8 6 /S 
Pi 5 Tle 4 Bindi eiG ee 
A he SAG | | 
“069 1634 1646 “eg a {aso 
ot S lew2Oule 46 | 
ete iwlc fe 2° sr assess = 


mm a ee ee He es es SS we tee ce PS es ES GP ES CS GS tee ee Se a ies ee ct Ee ee cee ES Ome ete Ome me me Ces oe ee eee whe tee me Oy ee ee ee 


ae a ama a a mh pe i iG QR pS ype: EE Se 





fe ae ete i + Se a a ee YR a en — - — 


ee = - —— te ee BS tee SS DS eeteny ny SE Ey es ee ey Sm ey ~ — eee re ee ee 











PROGQAM TO DETERMINE OOVARIANOB 


PUNOTION AND FOYER SPROTAUM 






T 


KOV = NO. SAMPLES 
[Le #8 








READ ENF(1), TL9S(1), TLOO(T), 
TU90(I), TU95(T) 


o> 
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OHG © (AK - EDF(J))/ 
(SDF (J 10 = EDF (J)) 
















L = TL9O(J) 


T9O 
(TL9O(J+1) - TL9O(J)) * CHG 






T95L = TL95(J) 
(TL95(J +1) = TLO5(J)) * OHG 


T95U = TU95(J) 
(TU95(J+1) - TU95(J)) * CHG 


T9OU = Waal 
(TUu90(J+1) - TU90(J)) * OHO 









Oo ® 





LOT 





108 






81 = 0.0 





$1 © 91+ O(K + 2) 
* OOS(tr* FK * FU/PH) 






J*#1)* 2.0 * 91 + O(1) 
+1) * cos(wv® Fu)/2 ° | 


+04 











U9OL = TOOL * U(T) 






V95L = TOSL * U(T) 






V9OU = TIOU * U(T) 








U95U = TI5U * UCI) 









WRITE! LAG, C(I-1) 
L(I), U9OL, U(I), UIOU, 
U95L, UCI), U9S5U 







OALL EXIT 





1G9 





LG 


"QOVAR" PLOTTING ROUTINE 


seal 


ISTAR = -15702990896 


IMINUS = 29809118256 ' 
IPLUS = 17997958192 
IBLANK = -17997953192 





Il = 27461634601 






“WRITI: 
HEADII'Gs 


2&0 
aye 


JOwi(1) = - } 
Jow2(I) = IMINUS | 
2 


0) - = of 


— 


WRITS 1 
DASH LINES 


290 
I = 2,52 


| 
| 


t 
ee 


JOW1(I) = IBLANK 


JOw2(I) = IKLANK 


290 
Jow1(1) = IY 
JOw1(53) = II . ; 
Jow2(1} = II 


JOw2(53) = II 





_ 8 
Ket 
° 













Ko UAE 


< 


oe >) 


SPANA = ALOGIO(UMAX * T9IOU) 
~ALOG1O(UMIN * T9OL) 







SFANB = ALOSIOCUMAX * T95U) 
-ALODIOC(UMIN * T9I5L) 


ZA & ALOGIO(UMIN * T9Ob) 


28 = ALOGIO(UMIN * TO5L) 


TA) = ((ALOG1O(U(T) * T9OL) 
~ZA)/3PANA) * 50.42. 







TA2 = ((ALO310(U(1))-za 
/3PANA * 50.+ 2. 


TA3} = ((ALOGIO(U(I) * T9OU) . 
~ZA)/SPAXA) * 50.4 2. 


JOWL( IAL) = IMINUS 
JOW(IA2) = ISTAR 


JOWICIA$) = IPLUS 


“((aLoa10(U(1) * T95L) 
-3B)/3PANB) * 50. 2. 


Ih) = 





182 = ((ALO910(U(1))-z8) 
/SPFAKB) ® SO. 2. 


IBS = ((ALOGIO(U(I) * T9SU) 
; -78)/SPANB) * £0. 2. 


JOd(THL) = IMINUS 
JEdz(1H7) = ISPAR 


JOUCTNS) = 


d 
(*) 


IPLUS 



















pelea 





JOVIC IAL) = 
JOWL(IA2) = 


JOV1(1A3) = 


JOVE(IK1) = 
JOW2( IR?) = 


JOW2(1B3) = 





IsLA’ K 


It LAK 


IBLAKK 


IFLADK 
I*LARK 


I4LAKK 





ie 
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APPENDIX D 


IRCA Me cays) 





sample X Three Millimeter Beads 
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Sample Y Two Size (4mm. and 5mm.) Beads 


* 





eis 
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APPENDIX E RESULTS OF AUTOCOVARIANCE PROGRAM i 


POWER SPECTRUM ESTIMATION 
———— FOR OLSCRETE TIME SERTES 
USING TUKEY-HANNEING WEIGHTING FACTORS 


TWO SIZE BEAOS ‘.40 ROWS MIXEO LAG = 40 
(1) 
SMOOTHED SPECTRUM WITH INDICATEO CONFIDENCE LIMITS 
90 PCT. 90 PCT. 95 PCT. 95 PCT. 
COVARIANCE RAW LOWER UPPER LOWER UPPER 
LAG SPECTRUM SPECTRUM LIMIT AVERAGE LIMtT LIMIT AVERAGE LIMIT 
0 0.9409 0.701 0.6647 0.7353 0.7911 0.6367 0.7353 0.7970 
l 0.6656 0.7684 0.6291 0.4959 0.7488 0.6027 0.6959 0.7544 
2 0.4372 0.5449 0.5591 0.6185 0.6655 0.5356 0.6185 0.6705 
3 0.2795 0.6159 0.5069 0.5607 0.6033 0.4856 0.5607 ® 0.6078 
4 “061668 0.4663 0.4365 0.4829 0.5196 0.4182 0.4829 0.5234 
5 0.0853 0.3831 0.3751 0.4149 0.4464 0.3593 0.4149 0.4497 
6 0.0401 0.4270 80.3815 0.4220 0.4541 0.3654 0.4220 0.4574 
u 0.0053 0.4507 | 0.35f0 0.3960 0.4261 0.3430 0.3960 0.4293 
a  -0.0057 "0.2556 0.2522 0.2790 0.3002 0-2416 0.2790 0.3024 
9 0.0084 0.1539 __ 0.1543 0.1707 0.1837 0.1478 0.1707 0.1851 
10 0.0375 0.1195 0.1187 0.1313 0.1413 0.1137 0.1313 021424 
ot 0.0604 061325 0.1165 9.1289 0.1387 0.1116 0.1289 0.1397 
12 0.057? 0.1310 0.1081 0.1195 0.1286 0.1035 0.1195 0.1296 
13 0.0475 0.0836 0.0951 0.1052 0.1132 0.0911 0.1052 0.1140 
14 0.0160 0.1227 0.0999 0.1105 0.1189 0.0957 0.1105 0.1198 
15 -0.0068 (041130 _ 0.0919 0.1017 0.1094 0.0880 0.1017 0.1102 
16 -0.0170 0.0580 0.0652 0.0721 0.0776 0.0624 0.0721 0.0782 
17 _-0.0260 — el poe Ge 0.0531 0.0588 0.0633 Meas 0.0509 0.0588 0.0637 
18 -0.0325 0.0563 0.0538 0.0595 “0.0640 °° &#420.0515 0.0595 0.0645 
19 -0.0214 | 0.0620 ~~ 0.0515 0.0569 0.0613 0.0493 0.0569 0.0617 
20 -0.0027 0.0455 0.0431 0.0477 0.0514 0.0413 0.0477 0.0517 
21 0.0354 __ 0.0380 0.0388 0.0430 0.0462 0.0372 0.0430 0.0466 
22 0.0533 0.0504 0.0418 0.0463 0.0498 0.0401 0.0463 0.0502 
23 0.0365 040464 | 0.0400 0.0442 0.0476 0.0383 0.0442 0.0479 
24 0.0242 0.0337 0.0361 0.0399 0.0430 0.0346 0.0399 ‘0.0433 
25 0.0232 0.0458 0.0366 0.0405 0.0436 0.0351 0.0405 0.0439 
26 0.0122 ~~ 0.0365 &# 0.0335 0.0370 0.0398 0.0321 . 0.0370 0.0401 
27 0.0137 0.0292 0.0268 0.0296 | 0.0318 0.0256 0.0296 0.0321 
28 0.0032 0.0235 0.0250 0.0276 °&#40.0297 0.0239 0.0276 0.0299 
29 _--0.0079 0 0343 0.0284 0.0314 0.0338 0.0272 0.0314 0.0340 
30 -0.0253 0.0335 ~ 0.0275 — 0.0304 } }#£46.0327 ~~ 0.0263 0.0304 0.0329. 
me) 70.0262 0.0204 _— 0.0223 0.0247 0.0266 0.0214 0.0247 0.0268 
32 -0.0133 “0.0246 £0.020% 0.0229 0.0246 ~~ 0.0198 ~~. 0.0229 “0.0248 
33 060067 0.0221  _— 00-0201 0.0222 0.0239 0.0192 0.0222 0.0241 
34 -0.0127 0.0201 0.0189 0.0209 °° — 0.0224 ~~ 0.0181 ~~ £0.0209 ~~ 0.0226 ~~ 
|..35 == 060162 == 000212 ———i200———s«é0-w 0221 0.0237 0.0191 0.0221 0.0239 
36 -0.0273 0.0258 0.0215 ~ 0.0238 ~ 020256 | 0.0206 © 060.0236 0.0257 ~~ 
__37_ . 7060333. 000222, 0.0221 __ 0.0244 0.0263 0.0211 0.0244 0.0265 
36 -0.0330 0.0274 0.0229 | 0.0254 © 0.0273 ~~ 0660220 = ———iO 0254 0275 : 
39. = 060391 000245 02090232 249 0.0201 0.0232 0.0251 
40 ~0.0413 0.0163 0.0184 “0.0204 #0.0219 ~— 6.0176 © 0.0204 0-0221 ° °°} © 
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POWER SPECTRUM ESTIMATION 
FOR OILSCREIF TIME SERIES 
USING TUKEY-HANNING WEIGHTING FACTORS 


TwO SIZE READS 40 ROWS MIXED LAG = 40 











(dD 
SMOOTHED POWER SPECTRUM SMOOTHEO POWER SPECTRUM 
WETH 90 PCT. CONFIDENCE LIMITS WITH 95 PCT. CONFIDENCE LIMITS 
0 l 2 3 4 5 6 7 8 9 0 0 I 2 3 4 5 6 7 8 9 0 

@eeaeeooaov env ep 2ceeeeseaeasneevneeeaevnee eee eee eeae@eeaaeaevneeevneae 2@eaoeaeteae eeeetspoaeaevpatepeoanvneeeaeveedgoeeaeeeeaevaeveeaeeeeaeaeeeeeeaee ea eaeseee bene ead 

0 I, o* ef 0 I : e e4¢{f 
ii 7 ¢ o #¢ | 1 I e of | 

2 1 0 e+ I 2i ; . * I 

3 | J ae ee ee e oF I 3 I o e+ I 
4 I 2 ae I 4 | ; e «4 I 
5 1 et ak . oF I 5 1 = . 04 I 

6 I : - er 2a I 6 I / « ¢ I 
71 oS ee oe I 71 7 + I 

8 I — ‘ . 0 I 8 I Ota : 
9 I 2 a4 I 9 I RS I 
10 f : esd, er 2 I 101 e oF I 
11 1 : as Se tlk : 2. oF I 
le | ae eee = de oe) 2 ee 7 oh 2 oF ...;.. eee I 
13 I 2 &¢ [I 13! o oF I 
14 I , - 7 ee | — | ,, i 
isin 0ti“(<C<=C;<C=Ci‘<‘“‘i‘i Oo ere. o a a CT Toe e + a I 
16 1. A nee ‘ : I 16 I « &¢ I 
i er a rs + he . oF we . | 
oi tsi‘(ititwtitit acer re ae me | ake e + ee > __l 
19 f- oe ¢ an t 19 (oe 7 . oF ee oe 7 i< 
201 _ 0 et ts) 3 = re, NIT ee I 20 I . e &4 I 
21 1 7 86 EN SRE ae tC 21 toe . oe piece tes —e ae 
ia | | «Ste. 2 ee.) 2 eee I 22 1 e+ I 
23 1 <—st >. oe ieee ac L233 fT aes ere a Se . ne I 
24i. e a4 I 24 I o 28s I 
25 1 704 1 251 . + ee ia 
26 I o o¢ I 26 I e 04 I 
i. 64° 7 aks. | sere rc | Paes | | a . a ra Y ween) OE EET SOOT TS GT STEED GR [a 
28 I oe I 28 I. o o4 I 
Pome 5  @6 0° eo . | ek ee eens) ee nS | er iia: 29 1 oan ony we ema 8S TEETER GGT UC~C~<“‘_‘“‘ ‘ ; ; ORS*‘“‘COSC™S i= 
30 I 2 &¢ I 30 1 e oF i 
—ai.s6 © © ee  @aemeenmgs ts SN Nm OO ec oo em Ge geese 
32 1 .e¢ I 32 1. e#¢ i 
———1.64  # =& °° °° °° + ©) oe A i. arn {35 1. 8 es ee eee eee gic + 2) an io 
34 [.0¢ 1 34 Ie o¢ I 
——.,¢¢ (ae () CergemiEn pao i | ae rt Snir iG Gp ia 
36 t .04 I 36 I .e¢ I 
ee 7+ | seus ge 3 ee 
38 tf .« o I 38 f . e+ I 
eee ee ee TY 39 1 eet SS = eee Ss eee ig ee 1G 
OC) ee [ 40 1.64 I 

2 TIME 1.3569 MIN ° ; 
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TwO STZE BEAOS 40 ROWS MIXEO LAG = 40 





(2) 
COVARIANCE RAW 
LAG SPECTRUM SPECTRUM 
0 ©3606 9748 0.6457 
1 0.7058 0.5820 
2 0.4705 0.7177 
> 0.2788 0.3926 | 
4 0.1308 0.7075 
ae 0.0277 0.5545 _ 
6 ~0.0383 0.5369 
| ue -0.0681 0.4480 
8 -0.0766 0.2408 
a5 -0.0500 0.2495 
10 -0.0153 0.0941 
11 0.0219 0.1406 
12 0.0435 0.0937 
13 0.0376 0.1299 
14 0.0379 0.0826 
15 0.0214 0.0902 
16 0.0149 0.0668 
_17_ 060164 0 0501 
|; 18 0.0054 0.0548 
we 19 0.0020 . 0.0537 
\ 20 -0.0065 0.0356 
21 —s—- - 0.0105 0.0500 
22 -0.0158 0.0343 
|_23 “0.0268  —— 0.0448 
me -0.0227 0.0333 
B25 -0.0250 0.0309 
| 26 -0.0252 ~ 0.0289 
|..27 70-0305 _——s_— 00358 
; 28 -0.0466 0.0250 
| 29 ~~ -0.0500_ 0.0358 
~ 30 -0.0460 ~~ 0.0196 
ol 70-0388 0.0219 
32 -0.0385 0.0164 > 
33 70001462 0.0358 
34 0.0189 “0.0276 
35 Ci 0533 ———— 0 0246 
36 0.0777 0.0220 
37. 060805, 0.0232 
38 0.0726 O.0175 © 
a9 _ 0.0557 —- 00256 
40 0.0262 0.0171 
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POWER SPECTRUM ESTIMATION 
FOR DISCRETE TIME SERIES 
USING TUKEY-HANNING WEIGHTING FACTORS 


SMOOTHED SPECTRUM WITH INDICATEO CONFIDENCE-LiFflts 


90 PCT. 
LOWER 
LIMIT 

0.5549 

0.5712 

0.5446 

0.4995 

0.5338 

0.5319 

0.4693 

0.3783 


0.2665 


0.1885 
0.1307 
0.1060 
0.1035 
0.0986 
0.0871 
0.0745 
0.0619 
0.0501 
0.0483 
0.0447 
0.0395 
0.0384 


0.0369 


0.0355 
0.0322 
0.0280 
0.0281 
0.0283 


0.0275. 


0.0262 


0.0219 


0.0180 
0.0204. 
0.0261 
0.0261 
0.0223 


~ 0.0208 


0.0194 


0.0189 


0.0194 


060247 


ere 


AVERAGE 
0.6138 
0.6318 
0.4025 
0.5526 
0.5965 
0.5884 
0.5191 
0.4184 
0.2948 
0.2085 
0.1446 
O-L1L?72 
0.1145" 
0.1090 
0.0963 
0.0824 
0.0685 
0.0555 
0.0534 
0.0495 
0.0437 
0.0424 
0.0408 
0.0393 
0.0356 
0.0310 
0.0311 
0.0313 


0.9364 


0.0290 
0.0242 
0.0199 
0.0226 
0.0289 
0.0289 


0.0230 
0.0215 
0.0210 
0.0215 


— 0.0214 


Ee a eo ee ee ee 


<2 & ~~ te 


90 PCI. 95 PCT. 
UPPER LOWER 
LIMIT LEMIT 

0.6605 0.5316 

0.6798 0.5472 

0.6483 0.5217 

0.5946 0.4785 

0.6354 0.5114 

0.6331 0.5095 

0.5585 0.4495 

0.4502 0.3624 

0.3172 0.2553 

0.2243 0.1805 

0.1556 001252 

0.1262 0.1015 

0.1232 0.0991 

0.1173 0.0944 

0.1036 0.0834 

0.0887 0.0714 

0.0737 0.0593 

0.0597 0.0480 

0.057% 0.0462 

0.0532 0.0428 

0.0470 0.0378 

0.0457 0.0368 

0.0439 0.0354 

0.0423 0.0340 

0.0383 0.0308 

0.0333 0.0268 

0.0334 0.0269 

0.0337 0.0271 

0.0327 0.0263 

0.0312 0.0251 

020261 0.0210 

0.0215 0.0173 

0.0243 5 0.0196 

0.0311 0.0250 

0.0311 “0.0250 

0.0266 0.0214 

0.0247 ~— 0.0199 

0.0231 0.0186 

a ae 0.0181 ~ 

0.0231 0.0186 

«9.02360 0°— 62-0185 


~~ "0.0225 


et rs a 2 


+ a eee we a ee 


AVERAGE 
0.6138 
0.6318 
0.6025 
0.5526 
0.5905 
0.5884 
0.4184 
0.2948 
0.2085 
0.1446 
0.1172 
0.1145 
0.1090 
0.0963 
0.0824 
0.0685 
0.0555 
0.0534 
0.0495 
0.0437 
0.0424 
0.0408 
0.0393 
0.0356 
0.0310 
0.0311 
0.0313 
0.0304 
0.0290 
0.0242 
0.0199 
0.0226 
0.0289 
0.02869 
0.0247 
0.0230 
0.0215 


~ 0.20210 


0.0215 


620214 


119 


95 PCT. 
UPPER 
LIMIT 

0.6654 

0.6849 

0.6531 

0.5990 

0.6401 

0.6378 

0.5627 

0.4536 

0.3195 

0.2260 

0.1567 

0.1271 

0.1241 

0.1182 

0.1044 

0.0894 

0.0742 

0.0601 

0.0579 

0.0536 

0.0474 

0.0460 

0.0443 

0.0426 

0.0386 

0.0336 

0.0337 

0.0340 

0.0329 

0.0315 

0.0262 

0.0216 

0.0245 

0.0313 


0.0313 — 


0.0268 


0.0249 — 


0.0233 
0.0227 
0.0233 
0.0232 
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POW!R SPECERUM ESTIMATION 
FCR MISCRETE LIME SERIES 
USING TUKFY-HANNING wWEFGHTING FACTORS 


TwO SIZE BEADS 40 ROWS MEXED LAG = 40 


(2) 
SMCOTHED POKER SPECTRUM : SMCOTHED PCWER SPECTRUM 
WITH 9C PCT. CONFINENCE LIMILIS WITH 95 PCI. CCNFIDENCE LIMITS 

0 l 2 3 4 5 6 fe R 9 0 0 1 2 3 4 5 6 7 8 9 0 
eeeoeeeoeeeeoaeteteoeaovseseeoeoeseeeesevanav,sv,evoeeee ses eeeeeoeececoeeceasvne7aeceeoaneae ee Beeeve @eeoeeoeoeeaneea ese oeoeeoaeoeaeoaeneeoeeoeae pees eoseeoeoeseoea ewe eoeaevn eva eceaseeanevn ev eae eee 
0 | eet [ 0 1 . * | 
' 1 oo +f lf « #6] 
2f ‘ e #¢ | 2 itl « 69 —f 
34 . ee f 31 o of | 
4 | o% + §f 4 1 o of fF 
Ss 4 o* + | 5 ] o @F | 
6 I 2 oF i 6 | 2 OF | 
ae | oF I ae | o 4 i 
Bf 084 I 8 | An I 
9 1 ot 4 | 9 { ry + l 
10 f 2 &4 l 10 f{ e + i 
ll 1 : 704 i it . 84 l 
12 1 ot i eee eee o #4 l 
13 { 2 + fF wt 2 oF H 
io { 2 OF t 4 o oF l 
15 { oe 4 i 19 | o #4 i 
16 1 oe 8+ i 16 { eo #4 [ 
i7‘l oe i? . * H 
18 1 of 4 H 18 { o 4 i 
19 1 2 #+ 1 19 | . #4 { 
20 { o H 20 { o a4 l 
H 764 1 211 2 8 l 
221 . 68 1 2215 2 88 | I 
23 { 0 #4 [ 23 ] eo e+ { 
i246 f o oF 1 24 { o &¢ . i 
Zo I oo 1 25 { e 4 i 
26 I o i 26 1 e + I 
ee? ft o @4 1 a7? f oe &4¢ | 
28 f o #4 i 281 o o+ I 
29 | a4 1 29 . 6 { 
30 T . i 30 f . #¢ i 
31 [eee fi 3l fees I 
32 1. oF 1 32 f 4 #¢ 1 
33 1 f+ i 33 {1 o + i 
34 1 284 1 34 | o #4 l 
35 — . o¢ i 35 § . «4 : H 
36 1. oF i 36 TT. «+ H 
37 1 #4 I 37 — e« # | 
3B 1. oF i 38 f. @F¢ ’ H 
39 | .e¢ 1 39 f. o¢ H 
40 I. #+¢ i 40 [. oF I 

2 TIME 123033 MIN 
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PCWFR SPECTHUM ESTIMATE IUN 
FOP OLSCREIF TiMe SERTES 
USING FLKEY-HANNENG WEIGHTING FACTORS 


 ; 
> 
o 


THREE SIZE PEACS 40 ROWS MIXED) LAG = 44) 
(3) 
SMGOTHEL SPECTRUM WITH INDICATED CONFIDENCE LIMITS 

90 PCI. 90 PCT. 95 PCI 95 PCI. 

COVAR LANCE RAw LCwtR UPPER LUWER UPPER 

SPFC FRUM SPECTRUM LOMit AVF RAGE LIM LiMit AVERAGE LiPis 

0 0.9035 0.45461 0.4371 0.4A35 0.5203 O.41A7 0.4835 0.5242 
l 0.6015 0.51t10 0.4342 0.4AhC3 0.5168 04-4159 0.48603 0.5206 
2 0.3676 0.4430 0.4367 Q.483] 0.5198 0.41863 C.4831 0.5237 
3 0e2Nt2 0.5353 0.4675 N.5172 0.5565 0.4479 0.5172 0.5606 
4 0.0977 0.5550 C.5005 9.5536 0.5957 0.4794 0.55346 0.6001 
5 0.0241 9.5691 1.4802 Ue5312 0.5715 0.4600 0.9312 0.575A 
6 -0.0311 0.4313 0.3973 0.4394 0.4728 0.3e&0€ 0.4394 0.476% 
T -0-0609 0. 3760 0.2918 0.3227 0.23473 0.2795 0.3227 0.3499 
8 “0.0816 0.2076 0.2131 0.7357 0.2536 0.2041 0.2357 90-2555 
9 -0.0742 0.7015 N.-1719 9.1902 22045 0.1647 0.1902 0.2062 
-0.0649 90-1502 0.1417 0.156" 0.1687 0.1358 0.1568 0.1699 
-0.0295 0.1251 0.1219 0.1348 0.1451 0.1168 0.1348 0.1461 
“0.0045 O- 1388 0.1261 0.1395 02.1502 0.1208 0.1395 0.1513 
0.0155 0.1554 0.1219 }. 1338 0.1440 0.1159 0.1338 0.1451 

0.9291 0.0857 0.N929 0-102A 0.1106 0.089C 0.1024 Oo L1I4 

0.0215 0.0844 O.NTAG 0-945 0.0909 0.0732 0.0845 0.0916 

0.0159 0.0835 0.0747 N.-OR26 | 0.0889 0.0715 0.0826 0.0896 

0.0212 0.0791 0.0776 0.9803 - 0.0A64% 0.0696 0.08CG3 0.0871 

0.0311 0.0797 _ 90-9642 0.0732 0.0788 0.0634 0.0732 0.0794 

0.0222 0.0546 0.0545 0.9603 0.0649 0.0523 0.0603 0.065% 

0.9099 0.0526 0.0477 0.0528 9.0568 0.0457 0.0528 0.0572 
-O0.N0A7 0.0513 0.0465 9.9515 0.0554 0.0446 0.0515 0.0558 
-0.0122 0.0508 0.0440 0.0487 0.0524 0.0422 0.0487 0.0528 
-0.0044 0.0420 0.036 0.9427 0.0459 0.037C 0.0427 0.0463 
0.0009 0.0361 0.0350 0.9387 0.0416 0.0335 0.0387 0.0419 

0.0088 0.0406 0.0376 0.9361 0.0388 0.0312 0.0361 0.0391 

0.01C3 0.0271 0.0315 0-9034R 0.0375 0.0302 0.0348 0.0378 

0.0014 0.0446 0.0353 0.03990 0.0420 0.033e 0.0390 0.0423- 
-0.0139 0.0397 0.0354 0.0391 0.0421 0.0339 0.0391 020424 
-0.0078 0.9325 0.03909 0.9342 0.0368 0.0296 0.0342 0.0371 
-0.0062 0.0320 0.0280 0.9310 0.0333 0.0268 0.0310 0.0336 
-0.0009 0.0273 0.0256 0.0283 0.0304 0.0245 0.0283 0.0306 
0.0053 120264 0.0246 0.9272 0.0293 0.0236 0.0272 0.0295 
-0.0058 0.02A7T 0.0243 0.9269 0.0289 0.0233 0.0269 0.0291 
-0.0093 0.0237 0.0247 0.9273 0-N294 0.023¢€ 0.0273 0.0296 
-0.0027 0.0331 0.0273 0.03G2 0.0325 0.0262 0.0302 0.0327 
-0.0024 0.0310 0.0272 0.0301 0.0324 0.0261 0.0301 0.0326 
~0.0123 0.0253 0.0240 0.9266 0.0286 0.0230 0.0266 0.0288 
-0.0208A 0.0247 0.0228 0.0252 0.0271 0.02186 0.0252 0.0273 
-0.0261 0.0261 020228 0.0253 0.0272 0.0219 0.0253 0.027% 
-0.036¢7 0.0241 0.0227 0.0251 0.0270 0.0217 0.0251 0.0272 
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PCWFER SPECIRUY ESTIMATION 
FCR OISCRETL TIME SERTES 
USIAG TUKEY-HANNING WEIGHTING 


40 ROWS MIXEO “LAG = 40 


SMCOTHED PCWER SPECTRUM 


WITH 90 PCT. 


0 I 2 


CONFIDENCE LIMITS 
4 5 6 7 4 9 0 
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FACTORS 


SMCOTHED PCWER SPECTRUM 
WITH 95 PCT. CONFIOENCE LIMITS 
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40 


COVART ANCE 


SPECTRUM 


0.92AL 
0.6461 
0.4014 
0.2204 
0.0973 
0.0130 
=e0'3 SE 
-0.9863 
-0.ORIA 
“0.0773 
-0.0 341 
0.0101 
0.0305 
9.0496 
0.0374 
0.0146 
0.00C6 
“0.0117 
-0.0219 
-0.04C9 
-0.03€I 
-0.0300 
-0.0171 
-0.0093 
-0.0001 
0.0065 
0.0055 
0.0096 
0.0225 
0.0228 
0.0294 
0.0209 
0.0073 
-0.0038 
-0.0L11 
-0.0184 
-0.0232 
-0.0141 
0.0002 


40 80WS MEXED 


RAW 


SPFC TRUM 


0.4511 
0.4A76 
0.5434 
0.5783 
0.5311 
0.6677? 
0.4A39 
0.2087 
0-1RB4 
0.2216 
0.1562 
0.152A 
0.1055 
0.0ROA 
0.0783 
0.0895 
0.0756 
0.0693 
0.0450 
0.9550 
0.0418 
0.0398 
0.0273 
0.945A 
0.0399 
0.0402 
0.0333 
0.025R 
0.0241 
0.0278 
0.0226 
0.032% 
0.0227 
0.0172 
0.0214 
0.0263 
0.0280 
0.0232 
0.0204 
0.0287 


USING 


POWER SPECTRUM ESTIMATION 
FER OLSCRITF TIME SERIES 


SMCOTHED) SPECTRUM WITH 


90 PCT, 
LOWER 
LIMIT 

0.4243 

0.4451 

0.4A6% 

0.5042 

0.5312 

0.43450 

0.28 14 

0.2024 

0.182% 

O17 

0.1552 

0.12A27 

0.1005 

0.07A0 

0.0739 

0.0754 

0.0701 

0.05%84 

0-04A% 

0.0445 

0.04C} 

0.04536 

0.031? 

0.0359 

0.0375 

0.0347 

0.0390 

0.0747 

0.0230 

0.9231 

0.0738 

0.0249 

0.0214 

0.0177 

0.0195 

0.0231, 

0.0238 

0.0214 

0.0209 

0.0222 


LAG 


=z 4() 


AVERAGE 
1.469% 
0.4924 
0.9 3R2 


"0.5578 


U.5771 
2.5876 
0.4812 
J- 3179 
0.2239 
0.7018 
0.1970 
Qeol7i? 
O.1L41R 
O.lill 
0.9863 
O.ORLR 
0.983% 
0.9776 
O.N64R 
9.9536 
0.9492 
0.0446 
0.9372 
0.90350 
0.9397 
0.930% 
0.9332 
0.9273 
0.9255 
0.0256 
0.0263 
0.0275 
0.9237 
0.0196 
0.0216 
0.9255 
0.0264 
0.0237 
0.9232 
0.0246 


TUKEY-HFANNING WEETGHTING FACTORS 


90 PCT. 
UPPER 
LIMIT 

9.25059 

0.5298 

0.5791 

0.6002 

0.6209 

0.6323 

0.5178 

0.3421 

0.2409 

0.2171 

0.2119 

0.1848 

0.1526 

0.1196 

9.0929 

0.ORRO 

0.0897 

).0R35 

0.0697 

0.0577? 

9.0529 

0.04RK0 

0.0400 

0.0377 

0.0427 

0.0446 

0.0413 

0.0357 

0.0294 

0.0274 

0.0275 

0.0283 

0.0296 

0.0255 

0.0211 


0.0232 — 


0.0275 
0.0284 
0.0255 
0.0249 
0.0264 


95 PCT. 
LOWER 
LIMEt 

0.4264 

0.466C 

0.4997 

0.5085 

0.4168 

0.2753 

0.1939 

0.1748 

0.17C6 

0.1487 

0.1228 

0.0962 

0.0748 

0.0708 

0.0722 

0.0€72 

0.0561 

0.0464 

0.0426 

0.038€ 

0.0322 

0.0303 

0.0344 

0.0359 

0.0333 

0.0287 

0.023€ 

0.0221 

0.0222 

0.0228 

0.0238 

0.0205 

0.0170 

0.0187 

0.0221 

0.0228 

0.0205 

0.0201 

0.0213 


INOICATEO CONFIOENCE LIPITS 


AVERAGE 
0.4694 
024924 
0.53e2 
0.557R 
0.5771 
2.5876 
O0.4R12 
9.3179 
9.2239 
0.2018 
0.1970 
O.1717 
O.141A 
O.1i11 
0.0863 
9.0834 
0.0648 
0.0536 
0.0492 
0.0446 
0.0372 
0.0350 
0.0397 
0.0415 
0.0384 
0.0332 
0.0273 
0.0255 
0.0256 
0.0263 
0.0275 
0.0237 
0.0196 
0.0216 
0.0255 
0.0264 
0.0237 
0.0232 
0.0246 


lize 


25 PC. 
UPPLR 
LIYMit 

0.5088 

0.5338 

0.5834 

0.6046 

0.6255 

0.6370 

0.5217 


‘0.3446 


0.242? 
0.2188 
Os 2135 
0.1861 
0.1537 
0.1205. 
0.0936 
0.0887 
0.0904 
0.0841 
0.0703 
0.05861 
0.0533 
0.0463 
0.0403 
0.0380 
0.0430 
0.0449 


0.0416 


0.0360 
0.0296 
0.0276 
0.0277 
0.0285 
0.0298 
0.0257 
0.0213 
0.0234 
0.0277 
0.0286 
0.0257 
0.0251 
0.0266 
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POWFR SPFCTRUM ESTIMATION 
FOR OISCRETE TIME SERIES 
USING TUKEY-HANNING WEIGHTING FACTORS 


THREF SIZE READS 40 ROWS MEXEO LAG = 49 
(4) 
SMOOTHEN POWER SPECTRIIM SMOOTHED PCWER SPECTRUM 
WITH 90 PCF. CONFIOENCE LIPIES WITH 95 PCT. CONFIDENCE LIMITS 
0 l 2 3 4 4 6 7 f 9 0 0 I 2 3 4 5 6 7 8 9 0 
e@eee¢eaecas eee eaeeCeaeaeaes eevee sespeeeCeae¢oasesevoseeaspeeene eee eseaagea ees cpeaeepesep eevee sseeoese eee ec oe se eee svpe esesaeseosvevoeeeeean eee see ees eaeeaevea 
0 | oat 1 0 | ° + 1 
| o o¢ I | ae | . ¢ I 
2! Pa I 2 i ° + I 
31 eo &¢ | 3 1 o + | 
41 es | 4 | . ee | 
5 1 2% +] 5 I . of 
6 1 o #4 I & I eo oF I 
7 1 . oF I ‘ame | ~ e4 I 
8 | . 86 I Se | « °° I 
9 1 8 4¢ I 9 1 « ot I 
10 1 oe6 I 101 I 
11 f oe¢ I tii! 1 
12 1 o o¢ I t2 1 I 
13 1 2 84 I 13 1 e I 
14 1 . &¢ I 14 ft o oF I 
15 I . oF I 15 1 o o¢ I 
16 1 0 0+ t 61 o o¢ I 
17 | o%4 I 17 4 o #4 I 
is I . oF I 186 I o o¢ l 
19 | 84 6 I 19 I ry + { 
20 |! « oF I 20 1 o oF I 
21 1 8 6 I 2iit a 4 06 1 
22 1 . 06 boee2.3 . oF I 
23 1 204 1 231 e¢ 1 
24 1 2 04 1 241 . oF I 
275 1 oe 1 25 1 2 06 1 
26 | oe I 261 « oF t 
27 1. Pa 1 27i . oF I 
28 1 7 oF I 28 1 7 oo I 
?9 «1 eo 8¢ I 29 1 o e¢ I 
30 (I e 96 I 30 I oe I 
31 1 204 | ee 2 ae | . #4 I 
32 I o o¢ I 32 1 « o¢ 1 
33 1 . OF I 33 1 « 84 I 
34 1.04 1 34 Te o¢ I 
35 — 04 1 35 fo #4 I 
36 I « OF I 361 . oF I 
37 1 084 I 37 I e o¢ I 
38 I .« oF 1 38 1 o oF 1 
39 I « 84 t 39 1 o 84 1 
40 1 oo 1 40 ft o 04 t 
2 «TIME 1.3056 MIN 
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POWER SPECTRUM ESTIMATION 
FOR OLSCRETE TIME SERIES 
USING TUKFY-<HANNING WEIGHTING FACTORS 


ome 


3mm BEAOS 40 RCWS MIXEO LAG > 40 
(s) 

_ SMCOTHED SPECTRUM WITH INOICATEO CUNFIOENCE LIMITS 
ie 90 PCI. 90 PCT. 95 PCT. 95 PCI. 
COVARLANCE RAW LOWER UPPER LOWER UPPER 
LAG SPEC TRUM SPECTRUM LIMLE AVERAGE LIMIf LIMit AVERAGE LIMIT 
0 0.8960 0.4985 0.4128 0.4567 0.4914 0.3555 0.4567 0.4950 
l 0.4939 0.4149 0.374h 0.4368 0.4700 0.3782 0.4368 0.4735 
2 0.2325 0.4190 0.3541 0.3917 0.4215 0.3392 0.3917 0.4246 
3 0.0761 0.3141 0.3093 0.34272 0.3682 0.2963 0.3422 0.3709 
4 0.0184 0.3215 0.3090 0.3418 0.3678 0.296C -0.3418 0.3705 
5 '0.0049 0.4102 0.3251 0.3596 0.3869 0.3114 0.3596 0.3898 
6 0.0114 0.2966 0.251 0.3154 0.3394 0.2731 0.3154 0.3419 
7 0.0276 0.2584 0.2473 0.2736 0.2944 0.2369 0.2736 0.2965 
8, 0.0287 0.2810 0.2475 0.2738 0.2946 0.2371 0.2738 0.2968 
9 0.0223 0.2749 0.2460 0.2721 0.2928 0.2357 0.2721 0.2950 
10 0.0125 0.2578 0.2307 0.755? 0.2746 0.2210 0.2552 0.2767 
11 0.0036 0.2305 0.2071 0.2291 0.2465 0.1984 0.2291 0.2483 
12 -0.0015 0.1976 0.1716 0.1898 0.2043 0.1644 0.1898 0.2058 
13 0.01C9 0.1337 0.1387 0.1534 0.1650 0.1328 0.1534 0.1663 
14 0.0237 0.1485 0.1237 0.1368 0.1472 60.1185 0.1368 0.1483 
15 0.0353 0.1166 0.1110 0.1228 0.1321 0.1C63 0.1228 0.1331 
16 0.0293 0.1093: 0.1006 U.1143 0.1198 0.0964 O.1113 0.1207 
17 0.0241 0.1101 0.0947 0.1048 07.1127 0.0907 0.1048. 0.1136 
18 0.0169 0.0896 0.0811 0.0897 0.0965 0.0776 0.0897 0.0972 
19 0.0054 0.0694 0.0704 0.0778 0.0838 0.0674 0.0778 0.0844 
20 0.0040 0.0830 0.0682 0.0755 0.0812 0.0654 0.0755 0.0818 
21 -0.0150 0.9665 0.0646 0.0714 0.0769 0.0619 0.0714 0.0774 
22 -0.0151 0.0697 0.0595 0.0658 0.0708 0.0570 0.0658 0.0714 
23 -0.0165 0.0574 0.0550 0.9608 0.0655 0.0527 0.0608 0.0660 
24 -0.0238 0.0589 0.0519 0.0574 0.0618 0.04S7 0.0574 0.0623 
25 -0.0298 0.0547 0.0475 0.0525 0.0565 0.0455 0.0525 0.0569 
26 -0.0144 0.0418 0.0416 0.0460 0.0495 0.0399 0.0460 0.0499 
27 -0.0053 0.0459 0.0415 0.0460 0.0494 0.0398 0.0460 0.0498 
28 0.0126 0.0502 0.0424 0.0469 0.0504 0.04Gé 0.0469 0.0508 
29 0.0179 0.0411 0.0379 0.0419 0.0451 0.0363 0.0419 0.0454 
30 0.0056 0.0352 0.0369 0.0408 "0.0439 0.0354 0.0408 0.0443 
31 0.0147 0.0519 0.0402 0.0445 0.0478 0.0385 0.0445 0.0482 
32 0.0188 0.0388 0.0397 90-0439 0.0472 0.038C 0.0439 02.0476 
33 0.0254 0.0460 0.0379 0.0420 0.0451 0.0363 0.0420 0.0455 
34 0.0202 0.0371 0.0357 0.0395 0.0425 0.0342 0.0395 0.0429 
35 0.0142 0.038! 0.0359 0.0397 0.0427 0.0344 0.0397 0.0430 
36 0.0119 0.0456 0.0397 0.0439 0.0472 0.0380 0.0439 0.0476 
37 0.0025 0.0462 0.0386 0.0427 0.0459 0.0370 0.0427 0.0463 
38 -0.0031 0.0327 0.0348 0.0385 0.0414 0.0333 0.0385 0.0417 
39 0.0128 0.0421 0.0348 0.0385 0.0414 0.0333 0.0385 0.0417 
40 0.0093 0.0370 0.0358 0.0396 0.0426 0.0343 0.0396 0.0429 
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POWFR SPECTRUM ESTIMATION _ 
FOR DISCRETE TIME SERTES a es 
USING TUKEY-HANNING WEIGHTING FACTORS 


34m BEADS 40 ROWS MIXED LAG = 40 
(6) 
SMCOTHEO POWER SPECTRUM SMCOTHEO POWER SPECTRUM 
WITH 90 PCT. CONFIOENCE LIMITS WITH 95 PCT. CONFIDENCE LIMITS 
O'. i r 4 3 4 5 6 7 8 9 (9) 0, 1 2 3 % 5 6 7 8 9 0 
eeeoecoeseoevoeeneooeae ep eeceoseesenoeseeoesepesseoevseeeeceoeceese eeeoaeaeeeeceaee eeceeeseoneeeoaosetP ose eee aoneneoaeegeoaesseosv boos teeaoavoteeseaeeves eee 
Oo 1  & of Oo 1 e eel 
1 iI o & + | 11 ry e¢ { 
21 - & @ 1 2! . 8 1 
ae | 2 24 I 3 1 . oF I 
& I . of I & | ee I 
5 I 78 ¢@ 1 5 I o oF 1 
6 I 7 © 4 1 6 1 « oF { 
7 1 7 © + i rit « oF I 
6 I 7 & ¢ 1 8 I . 1 
9 |! e & 4 if 9 1 ry e¢ I 
10 1 . & + 1 tof « %4 1 
11 1 7 # ¢ 1 1 « oF 1 
12 I 9 4 i tet « oe i 
13 1 oe + 1 31 . oF { 
14 I « & ¢ 1 14 . 8 I 
15 1 8 ¢ 1 15 . oF 1 
16 I 28 ¢ 1 61 « 8 1 
17 { « * ¢ [ 17 { e o¢ I 
18 { « & 6 : “i J ie f ° oe { 
veld f o #4 | 19 1  * 4 1 
70 1 o Fe + 1 20 { . 8 { 
211 7 8 + 1 21 7 «(86 : 1 
22 1 o 8 i 22a . 8 + { 
73 1 o & ¢ J 23 { e e¢ 1 
24 1 oe %¢ 1 24 1 . 86 1 
25 1 - & 4 1 25 { 2 8 1 
26 { 7 1 261 . oF I 
~ 27 | « 8 2791 . oF I 
26 I  & + { 28 I . oF I 
Vine 9 Oe avai Se i. 1 29 1. e4 { 
30 { . a6 I 30 Te @ I 
YIr.e¢ a - t  3t ft 2. «F al 
32 = « 8 { 321 . #¢ { 
33 1. 2 + ca . { 33 t . o¢ : ji 1 
34 ft. o¢ 1 34 Ie e¢ { 
35 I. 2 @ 1° 35°Ts -e¢€ a at 
36. 1 es OF 1 36% . #e¢ i 
ae ee es UC” a ria oa. 1 3% ts oe eee Se a an sae ae et 
38 [ee + | ef ea { 
39 1.9 @ ‘<i 139%. 8 @ i 7 me" Oh 
oe a 1 40 1. 0 ' 
i EE GR Jel 
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PCWER SPECTRUM ESTIMATION 
FOR OISCRETF TIME SERIES 
USING TUKEY-HANNING WEIGHTING FACTORS 


3MM BFADS 40 RCWS MIXEO LAG = 40 
(6) 
SMGOTHED SPECTRUM WITH INDICATED CONFIDENCE LIMITS 

90 PCT. 90 PCT. 95 PCT. 95 PCT. 

COVAR TANCE RAW LOWER UPPER LOWER UPPER 

LAG SPEC TRUM SPECTRUM . LIM(ET AVERAGE LIMIT LIMLy AVERAGE LIM 
0 0.937f/ 0.3632 0.336? 0.3725 0.4008 0.3225 0.3725 0.4037 
1 0.5903 O.3AL? 0.3256 0.3602 0.3876 0.3119 0.3602 0.3904 
2 0.2076 0.3141 0.3395 0.3755 0.4041 0.3252 0.3755 0.4071 
3 0.01€0 0.4923 0.36597 0.4045 0.4352 0.3503 0.4045 0.4385 
4 -0.0592 0.3193 0.3115 0.3446 0.3708 0.2984 0.3446 0.3735 
5 -0.04C6 0.2473 0.252n 0.2788 0.2999 0.2414 0.2788 0.3022 
6 0.0211 ~O.3011 0.253? 0.2807 0.3020 0.2431 0.28C? 0.3042 
7 0.07S7 0.2732 90-2560 0.2832 0.304? 0.2453 0.2832 0.3070 
M 0.0aC9 0.2854 0.2612 0.28R9 0.3108 0.2502 0.2889 0.3132 
9 0.0265 0.3115 0.2886 0. 3193 0.3435 0.2765 0.3193 0.3461 
10 -0.0284 0. 3686 0.3366 0.3723 0.4006 0.3224 O. 3124 0.4036 
lt -0.0670 0.4406 0. 3343 0.3698 0.3979 0.3202 0.3698 0.4009 
12 -0.0634 0.7293 0.2334 0.7582 0.2778 0.2236 0.2582 0.2799 
13 -0.0356 0.1336 0.1479 0.1636 0.1761 0.1417 0.1636 0.1774 
14 -0.0014 0.1580 0.1246 0.1378 0.1483 0.119% 0.1378 0.1494 
15 0.0178 0.1016 0.1045 0.1156 0.1244 0.1col 0.t156 0.1253 
16 0.0163 0.1013 0.0976 0.1025 0.1103 0.0887? 0.1025 O.-liti 
7? -0.0261 0.1057 0.0835 9.0973 0.0994 0.0800 0.0923 0.1001 
18 -0.0426 0.0566 0.0692 0.0765 0.0824 0.0663 0.0765 0.0830 
19 -0.0340 0.0872 0.0719 0.0817 0.0879 0.07C8 0.0817 0.0886 
20 -0.0078 0.0959 0.0744 0.0823 0.0886 0.0713 0.0823 0.0892 
21 0.0252 0.0503 0.0669 0.0731 0.0786 0.0633 0.0731 0.0792 
22 0.0457 0.0959 0.0636 0.0703 0.075? 0.0609 0.07C63 0.0762 
23 0.0485 0.0394 0.0544 0.0602 0.0648 0.0521 0.0602 0.0653 
24 0.0337 0.0662 0.0529 0.0585 0.0629 0.0506 0.0585 0.063% 
25 0-016” 9.0620 0.0552 0.0610 0.0657 0.0529 0.0610 0.0662 
26 0.0096 0.0538 0.0481 0.0532 0.0573 0.0461 0.0532 0.0577 
27 0.0069 0.0432 0.0431 0.0477 0.0513 0.0413 0.0477 0.0517 
28 0.0445 0.0504 0.0426 9.0471 0.050? 0.0408 0.0471 0.0511 
29 0.0595 0.0446 0.0429 0.0474 0.0510 0.0411 0.0474 0.0514 
40 0.016&4 0.0502 0.045R 0.0506 0.0545 0.04386 0.0506 0.0549 
31 -0.0203 0.0575 0.0465 0.0514 0.0553 0.0445 0.0514 0.0557 
32 -0.0226 0.0403 0.0415 9.0459 0.0494 0.0397 0.0459 0.0497 
33 -0.0139 0.9453 0.0409 0.0452 0.0487 0.0392 020452 0.0490 
4 -0.02C64 0.9500 0.0425 0.0470 0.0506 0.0407 0.0470 0.0510 
35 -0.0295 0.0429 0.0411 0.0455 0.0489 020394 0.0455 0.0493 
36 -0.0031 0.0462 0.0386 0.0427 0.0459 0.0369 0.0427 0.0462 
37 -0.01 34% 0.0354 0.0367? 0.0405 02-0436 0.0351 0.0405 0.0440 
38 -0.0295 0.0452 0.0369 0.0408 0.0439 0.0354 0.0408 0.0443 
39 -0.0343 0.0375 0.0398 0.0441 0.0474 0.0382 0.0441 0.0478 
40 -0.0159 0.0560 0.0423 0.0468 0.0303 0.0405 0.0468 0.0507 
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PCWER SPECTRUM ESTIMATION 
FOR DISCRETE*TIME SERTES 


USING TUKEY-HANNING WEIGHTING FACTORS 


a 


34M BEADS 


apt mw ee nee 


40 RCWS MIXED LAG = 40 
(6) 
SMCOTHED POWER SPECTRUM 
WITH 90 PCT. CONFIDENCE LIMITS 
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SMOOTHED PCWER SPECTRUM 
WITH 95 PCT. CONFIOENCE LIMIT 
1 2 3 4 - 6 Ts 


128 


8 9 


0 


@eeeoeeoeseosesoeoeeoe eevee ese eoeoeeeeoeaeceantoseeetsteoeteoeeseaoseenenesd 


e¢ 


e o¢ 
as 
. 06 
e a4 . 
e oe ¢ 
e+ 
e¢ 
e ¢ 
ee 
e & @ 
e ¢ 
oe 
. oe 
° e 4 


e e+ 
ny’ 
e a¢ 
P e* 
e + 
. of 
e @+ 
P + 
° e+ 
o o4 
.- of 
ry a+ 
Le 
| 
| 
I 
l 
I 
l 
| 
| 
J 
m 
o I 
| 
et 
| 
eee. 
| 
I 


Sees ne 


tao ape 


te 


= Tew eT Te ae ee ee es ee ee oe eee a ee oe 
- 


os 


























hes ree Tenge 





ill I tld | il il 


init 


